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ABSTRACT 


In  this  thesis,  the  Monterey-Miami  Parabolic  Equation  (MMPE)  model  is  used  to 
generate  predictions  from  numerical  analysis  of  the  reverberation  loss  structure  and  peak 
vertical  correlation  structure  generated  by  the  water/bottom  interface,  the  bottom/sub¬ 
bottom  interface,  and  the  bottom  volume  for  a  shallow  water  environment.  These 
predictions  are  then  compared  to  the  peak  vertical  correlation  analysis  of  recorded  data 
collected  in  an  actual  shallow  water  environment  similar  to  the  modeled  environment. 
This  experimental  data  was  recorded  by  a  32-element  vertical  line  array  (VLA)  that 
recorded  the  reverberant  return  generated  by  charges  detonated  over  the  continental  shelf 
in  the  East  China  Sea  as  part  of  ASIAEX.  A  comparison  is  made  between  predictions 
and  recorded  data  by  analyzing  trends  in  peak  vertical  correlation  with  decreasing 
bandwidth.  The  influences  of  interface  roughness,  bottom  volume  perturbations,  and 
water  volume  turbulence  on  peak  vertical  correlation  is  also  detennined. 
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I.  INTRODUCTION 


Amazing  advances  have  been  achieved  in  the  implementation  of  tactical  active 
sonar  systems  to  prosecute  the  underwater  environment  since  their  advent  over  eight 
decades  ago.  In  all  of  this  time,  however,  very  little  progress  has  been  made  in 
combating  the  adverse  effects  that  strong  reverberation  in  a  shallow  water  environment 
has  on  the  effectiveness  of  these  systems.  As  the  United  States  Navy  places  an  ever- 
increasing  effort  on  operating  in  the  littoral  regions  of  the  world,  finding  methods  to 
understand,  predict,  and  compensate  for  the  effects  of  reverberation  to  increase  the 
performance  of  underwater  detection  and  tracking  systems  in  the  shallow  water 
environment  has  become  more  important  than  ever.  To  realize  such  an  advance  requires 
a  better  understanding  of  the  primary  mechanisms  that  generate  reverberation  in  the 
shallow  water  environment,  the  impact  of  environmental  variability  on  acoustic 
propagation,  and  the  unique  factors  of  the  reverberation  returns  that  may  allow  signal 
discrimination  by  advanced  processing  techniques. 

To  study  the  effects  of  deep  ocean  acoustic  reverberation,  a  series  of  experiments 
was  conducted  near  the  Mid-Atlantic  Ridge  over  a  several  year  span  in  the  early- 1 990 ’s. 
This  program  was  known  as  the  Acoustic  Reverberation  Special  Research  Project 
(ARSRP)  and  was  sponsored  by  the  Office  of  Naval  Research  (ONR).  The  results  of  the 
data  analysis^'  [3]'  |4J  from  these  experiments  showed  that  the  general  structure  of  the 
measured  reverberation  often  agreed  with  the  predictions  that  were  based  on  propagation 
modeling,  particularly  when  high  reverberant  returns  were  observed  in  areas  where  the 
acoustic  energy  interacted  strongly  with  the  bottom.  This  conclusion  validated  the  use  of 
propagation  modeling  to  accurately  predict  the  observed  long  wavelength  reverberation 
statistics  of  the  two-way  transmission  loss. 

Recently,  another  such  study  shifted  the  focus  of  reverberation  research  to  the 
shallow  water  environment.  ONR  sponsored  this  reverberation  study  as  part  of  the  Asian 
Seas  International  Acoustic  Experiment  (ASIAEX),  a  scientific  collaboration  between  the 
United  States,  the  People’s  Republic  of  China,  and  Japan,  conducted  in  June  of  2001. 
This  experiment  was  conducted  over  the  shallow  continental  shelf  of  the  East  China  Sea. 
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In  this  region,  a  soft  bottom  layer  of  sand  or  mud  lies  about  100m  below  the  ocean 
surface,  with  a  harder  layer  of  sand,  gravel,  or  rocks  beneath  it.  Evidence  suggests  that 
that  the  transition  between  these  two  layers  fonns  a  second  distinct  interface  layer  that 
lies  approximately  10m  below  the  water/bottom  interface.  Recorded  data  from  the 
reverberation  generated  by  explosive  sources  and  recorded  by  a  32-element  vertical  line 
array  (VLA)  during  this  study  is  available  thanks  to  the  efforts  of  Ji-Xun  Zhou  of  the 
Georgia  Institute  of  Technology. 

In  recent  years,  continuing  thesis  work^’  l?l  has  been  conducted  using  the 
Monterey-Miami  Parabolic  Equation  (MMPE)  propagation  model1*1  to  predict  the 
structure  of  reverberation  in  a  shallow  water  environment.  Smith  and  Li  first  introduced 
the  use  of  this  method  in  2000. [21  Since  then,  the  model  has  been  updated  and  improved 
to  more  accurately  simulate  the  ocean  environment  and  the  effect  of  its  features  on 
acoustic  propagation.  Such  improvements  now  accurately  predict  the  effects  of  bottom 
interface  roughness,  bottom  volume  perturbations,  and  water  volume  turbulence. 

The  previous  thesis  work  involved  a  modeled  environment  and  geometry  similar 
to  that  observed  during  the  ASIAEX.  In  this  work,  a  16-element  VLA  was  modeled  to 
support  predictions  of  both  monostatic  and  bistatic  reverberation  in  a  shallow  water 
environment.  The  element  at  a  depth  of  48m  was  used  to  act  as  a  source  placed  in  the 
middle  of  the  water  column  and  all  of  the  elements  acted  as  receivers.  The  water/bottom 
interface  was  set  at  100m  with  a  bottom/sub-bottom  interface  placed  5m  below  it.  The 
modeled  reference  environment  included  bottom  and  sub-bottom  characteristics  based  on 
geo-acoustic  data  taken  from  the  East  China  Sea  provided  by  James  H.  Miller  and  Gopu 
Potty  of  the  University  of  Rhode  Island  and  a  sound  speed  profile  derived  from  averaged 
Current  Temperature  Depth  (CTD)  measurements  taken  from  the  R/V  Melville  provided 
by  Peter  Dahl,  Chief  Scientist  for  the  East  China  Sea  component  of  ASIAEX.  Within 
this  modeled  environment,  reverberation  analyses  were  conducted  for  continuous  wave 
(CW)  pulses  and  broadband  pulses  for  a  variety  of  different  environmental  variations.  To 
generate  these  broadband  pulses,  the  MMPE  Model  was  run  for  512  frequency  bins 
spanning  a  250  Hz  bandwidth  centered  on  a  frequency  of  250  Hz.  In  these  analyses, 
particular  attention  was  paid  to  reverberation  loss  generated  by  the  interfaces  and  the 
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bottom  volume  and  the  peak  vertical  correlation  of  these  reverberation  loss  structures. 
Some  of  the  significant  findings  from  these  previous  works  are  as  follows: 

•  CW  analysis  was  unable  to  capture  the  coherent  structure  of  reverberation 
from  the  bottom  volume  due  to  an  inability  to  properly  resolve  multi-path 
effects.  As  a  result,  peak  vertical  correlation  analysis  was  only  valid  for 
broadband  calculation. 

•  Peak  vertical  correlation  analysis  showed  that  reverberation  from  the 
volume  decorrelates  more  rapidly  across  the  vertical  array  than 
reverberation  from  the  water/bottom  interface  and  the  bottom/sub-bottom 
interface. 

•  Higher  reverberation  levels  are  observed  with  increased  roughness  of  the 
interface.  Peak  vertical  correlation  for  the  bottom/sub-bottom  interface 
displayed  a  clear  decreasing  trend  with  increased  roughness. 

•  Lower  reverberation  levels  are  observed  for  larger  sound  speed 
perturbations  in  the  bottom  volume.  Despite  the  presence  of  this  trend,  no 
significant  change  is  observed  for  the  peak  vertical  correlation  for  either 
interface  or  the  bottom  volume  with  larger  sound  speed  perturbations. 

•  The  influence  of  variations  in  the  sound  speed  profile  (SSP)  was 
insignificant  with  respect  to  both  reverberation  levels  and  the  peak  vertical 
correlation  for  either  interface  or  the  bottom  volume. 

The  objective  of  this  analysis  is  to  expand  on  the  previous  work  and  use  the 
methods  developed  in  them  to  generate  accurate  reverberation  predictions  for  the 
environment  observed  in  ASIAEX.  In  conducting  the  numerical  analysis,  the  source  will 
be  modified  from  previous  work  to  more  accurately  model  the  explosive  charges  used  in 
ASIAEX.  To  account  for  a  larger  portion  of  the  acoustic  energy  generated  by  a  wide 
band  explosive  source,  a  wider  frequency  bandwidth  of  500  Hz  centered  at  350  Hz  as 
calculated  over  1024  frequency  bins  will  be  used.  The  frequency  spectrum  generated  by 
a  Standard  Underwater  Source  (SUS)  charge,  similar  to  the  charges  used  in  ASIAEX, 
will  also  be  applied  to  the  numerical  analysis  instead  of  the  Hanning  window  frequency 
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spectrum  applied  in  the  previous  work.  This  will  be  accomplished  by  applying  a 
modeled  SUS  charge  spectrum  as  provided  by  David  Knobles  of  the  University  of  Texas. 

The  results  of  the  predictions  generated  by  the  numerical  analysis  will  then  be 
used  as  a  basis  for  comparison  with  the  recorded  data  from  ASIAEX.  Recent  research'131 
describes  a  method  in  which  the  spatial  coherence  of  a  reverberant  signal  can  yield 
information  regarding  the  nature  of  the  reverberation  -  specifically,  whether  this 
reverberation  is  generated  at  an  interface  or  within  the  bottom  volume.  This  method 
exploits  changes  in  spatial  coherence  with  changes  in  transmitted  pulse  length,  which 
implies  that  similar  changes  can  be  exploited  for  different  bandwidths  when  using  a 
broadband  source.  As  a  measure  of  spatial  coherence,  emphasis  will  be  placed  on  the 
peak  vertical  correlations  for  both  the  predicted  and  recorded  signals  and  how  these  peak 
vertical  correlations  are  effected  by  bandwidth  size. 

Exploitable  changes  in  spatial  coherence  are  those  that  display  some  type  of 
distinct  difference  between  the  interfaces  and  the  bottom  volume  in  the  behavior  of  their 
peak  vertical  correlation  structures  as  bandwidth  is  changed.  Variations  in  environmental 
factors  could  affect  the  coherence  of  the  acoustic  signal  through  both  direct  effects  on 
scatter  and  indirect  effects  along  the  propagation  path.  If  the  coherence  is  significantly 
affected,  a  change  may  then  arise  in  peak  vertical  correlation  structures,  which  could 
diminish  or  destroy  exploitable  changes  in  spatial  coherence.  To  detennine  the  influence 
of  direct  effects  on  scatter,  variations  in  roughness  for  both  the  water/bottom  interface 
and  the  bottom/sub-bottom  interface  and  variations  in  the  magnitude  of  bottom  volume 
sound  speed  perturbations  are  generated  and  the  resulting  reverberation  loss  and  peak 
vertical  correlations  are  determined.  To  examine  an  indirect  effect,  variations  in  water 
volume  sound  speed  turbulence  will  be  examined  in  the  same  fashion.  Similar  analysis 
for  the  direct  influences  of  interface  roughness  and  bottom  volume  sound  speed 
perturbations  has  been  done  in  previous  work,'-71  but  not  for  the  source  spectrum  that  will 
be  applied  here. 

The  remainder  of  this  thesis  is  organized  as  follows.  The  second  chapter 
introduces  reverberation  theory,  environmental  variation,  and  the  use  of  spatial  coherence 
for  separating  and  discriminating  bottom  and  surface  scatter.  Methods  for  the  numerical 
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implementation  of  this  theory  will  be  discussed  in  the  third  chapter.  The  results  of  both 
the  numerical  analysis  and  analysis  of  the  recorded  data  from  ASIAEX  and  their 
implications  are  presented  in  the  fourth  chapter.  The  final  chapter  summarized  this  work 
and  generates  a  few  suggestions  for  future  work. 
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II.  THEORY 


A.  REVERBERATION  THEORY 

For  underwater  acoustics  applications,  the  term  “reverberation”  pertains  to  the 
scattering  and  re-radiation  of  acoustic  energy  transmitted  from  a  source  by  the  multitude 
of  non-target-of-interest  items  in  the  ocean  environment.  These  items  include  waterborne 
elements  such  as  bubbles,  marine  life  (the  swim  bladders  of  fish  in  particular),  and  the 
ocean  surface.  Elements  of  the  ocean  bottom  such  as  the  sea  floor  interface  and 
discontinuities  in  the  bottom  volume  also  play  a  significant  factor.  In  the  shallow  water 
environment,  it  is  the  effects  of  this  later  group  that  is  of  particular  importance.  To  better 
understand  the  mechanisms  by  which  bottom  features  contribute  to  reverberation,  the 
effect  of  the  water/bottom  interface,  a  bottom/sub-bottom  interface,  and  inhomogeneities 
within  the  bottom  volume  are  examined. 

The  scattering  mechanisms  for  reverberation  for  both  a  monostatic  geometry,  in 
which  source  and  receiver  are  collocated,  and  a  bistatic  geometry,  where  source  and 
receiver  are  not  collocated,  are  identical.  The  following  development  describes  the 
monostatic  geometry  (source  and  receiver  may  be  separated  in  depth  only)  since  its 
numerical  implementation  is  somewhat  simpler. 

Mean  reverberation  pressure  level,  RPL,  can  be  defined  as15 


RPL  =  SL  +  DIt  +  DIr  + 10  log 


M, 

Rl 


-RL, 


'b,v  9 


(1) 


where  SL  is  the  source  level,  DIT  is  the  directivity  index  for  the  transmitter,  DI R  is  the 
directivity  index  for  the  receiver,  AAh  is  the  ensonified  area  in  the  horizontal  direction, 
R0  is  the  reference  distance,  and  RLb  v  is  the  reverberation  loss  per  unit  area  for  either  the 

bottom  interface  or  the  volume,  respectively.  Where  appropriate,  these  terms  have  units 
of  decibels  relative  to  1  (iPa  and  the  reference  distance  of  1  m.  At  this  portion  of  the 
theoretical  development,  reverberation  loss  per  unit  area,  RLh  v ,  (or  simply  reverberation 
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loss)  is  a  full-wave,  CW  quantity  and  does  not  explicitly  depend  on  the  angle  of 
interaction.  However,  because  the  scattering  model  is  based  upon  the  parabolic 
approximation  in  the  far-field,  a  near  grazing  angle  of  interaction  is  still  assumed. 

1.  Bottom  Interface  Scattering 

Reverberation  loss  per  unit  area  for  the  bottom  interface,  RLh ,  is  defined  as[5] 

M,=-201og(|ft  (2) 

where  p2_wayb  is  the  average  nonnalized,  two-way  pressure  field  (based  on  the  long 

wavelength  components  of  environmental  fluctuations)  from  the  transmitter  to  the 
scattering  patch  at  the  bottom  and  back  to  the  receiver  (to  be  defined  in  the  next  chapter) 
and  Sb  is  the  full-wave  scattering  strength  from  small-scale  interface  roughness.  The 

interface  scattering  strength,  Sb ,  is  defined  as[5] 

S„  =  101og  -iy  —  Uf(2i,)  ,  (3) 

32 n-  \  c0  ) 

where  k0  is  the  reference  wavenumber,  Ac  is  the  difference  in  sound  speed  between  the 

water  and  the  bottom  at  the  interface,  c0  is  the  reference  sound  speed,  and  W2\lkQ )  is 

the  two-dimensional  (2-D)  spectrum  of  the  interface  roughness  evaluated  at  the  Bragg 
wavenumber  for  monostatic  reverberation.  The  distinction  between  long  and  short 
wavelength  components  is  chosen  to  be  on  the  order  of  a  few  acoustic  wavelengths. 

2.  Volume  Scattering 

The  nature  of  the  volume  reverberation  loss  is  more  complex,  and  it  cannot  be 
treated  simply  in  terms  of  transmission  loss  adjusted  by  a  scattering  term.  It  can, 
however,  be  defined  as[5] 

RLy  =  -20 log  j  p2.wayy(z)dz  -Sv,  (4) 

jb+ni 
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where  p2_  (z)  is  the  normalized  two-way  pressure  field  within  the  volume  (to  be 

defined  in  the  next  chapter)  and  Sv  is  the  volume  scattering  strength.  Similar  to  the 
interface  scattering  strength,  Sh ,  the  volume  scattering  strength,  Sv ,  is  defined  as[5] 


Sv  =  10  log 


32  n 


2  <} 


(2  k0) 


(5) 


where  W2H]  (2  k0)  is  the  2-D  horizontal  spectrum  of  the  volume  fluctuations.  This  later 
factor  is  assumed  to  be  horizontally  isotropic  and  independent  of  depth  when  evaluated  at 
the  Bragg  wavenumber. 

It  is  important  to  note  that  all  terms  other  than  RLbv  in  Eq.  (1)  are  treated  as 
constants.  The  scattering  strengths,  Sb  v ,  are  also  assumed  to  be  generally  constant  over  a 

region  of  statistically  similar  perturbations.  Therefore,  the  structure  of  the  mean 
reverberation  pressure  level,  RPL,  is  dominated  by  the  propagation  defined  by  the  two- 
way  normalized  pressure,  p2_way  {h  v} . 

B.  INTERFACE  PERTURBATION 

For  acoustic  scattering  from  the  interface,  it  has  been  shown  that  reverberation 
loss  for  the  bottom  interface,  RLh ,  is  dependent  upon  the  bottom  scattering  strength,  Sh 

(Eq.  (2)).  This  scattering  strength  is,  in  turn,  dependent  on  the  two-dimensional  (2-D) 
spectrum  of  the  interface  roughness  evaluated  at  the  Bragg  wavenumber  for  monostatic 

reverberation,  W2\2k0)  (Eq.  (3)).  To  further  define  this  2-D  interface  roughness,  we 
assume  a  two-dimensional  (2-D)  interface  spectrum  of  the  form'-1  '' 

w!',(K)  = - — sr.  (0 

d+4 XV1 

where  kr  is  the  horizontal  spatial  wavenumber  vector  defined  by  kr  =  ^K2  +  L2  ,  K  and  L 
are  the  horizontal  wavenumbers  in  the  x-  and  y-directions,  respectively,  //  is  a 

nonnalization  factor,  Lcorr  is  a  correlation  length  scale,  and  (3  is  the  spectral  exponent.  If 
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this  2-D  interface  spectrum,  W2]  (kr) ,  is  also  assumed  to  be  isotropic,  then  the 

nonnalization  factor,  /u,  can  be  defined  in  terms  of  the  root-mean-square  (nns)  roughness, 
<7,  such  that 


where 


2k\  W(2I\kr)krdkr=CJ2 , 

o 


M  = 


n 


P 


^-1 

v  2 


_2  T  2 
0  L 

corr 


(V) 


(8) 


The  scattering  amplitude  can  then  be  determined  by  evaluating  W2(1)  (kr )  at  kr=2k0  since 

it  is  the  result  of  Bragg  scatter  as  evaluated  along  the  line  of  propagation  for  monostatic 
reverberation. 


To  develop  an  expression  for  the  long- wavelength  interface  roughness,  the  one¬ 
dimensional  (1-D)  spectrum  along  the  x-axis  is  also  required.  This  can  be  accomplished 
by  taking  the  1-D  transform  of  the  2-D  interface  spectrum,  W2!)  ( K,L ),  along  a  slice  at 
y=0,  which  yields 

WlI]{K )  =  j  W2\k  ,  L)dL  =  Y(J2Lcon.(\  +  L2corrK V  "  (9) 

where 


r= 


P-1 

r\ 

1 1 

2  J 

UJ 

2 


jtr 


(P) 


(10) 


C.  VOLUME  PERTURBATION 

1.  Sound  Speed  Fluctuations 

Similar  to  acoustic  scattering  from  the  bottom  interface,  reverberation  loss  due  to 
the  bottom  volume,  RLV ,  is  dependent  upon  the  volume  scattering  strength,  Sv  (Eq.  (4)). 
In  turn,  the  volume  scattering  strength  is  dependent  on  the  2-D  horizontal  spectrum  of  the 

volume  fluctuations,  W2H\2k0 )  (Eq.  (5)).  To  attain  this  2-D  spectrum,  the  volume  sound 
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speed  perturbations  are  first  modeled  by  using  a  three-dimensional  (3-D)  volume 
spectrum  given  by[12] 

W\v\k,L,M)  =  ^^-[a1(K1 +I})  +  M2J^~1 ,  (11) 

2tc 

where  B  is  the  spectral  strength  constant,  /?  is  the  spectral  exponent,  A  =  —  =  —  is  the 

horizontal-to-vertical  aspect  ratio  used  to  account  for  the  anisotropy  of  fluctuations  in  the 
sediment,  K  and  L  are  the  horizontal  wavenumbers  in  the  x-  and  y-directions, 
respectively,  and  M  is  the  vertical  wavenumber. 

From  this  3-D  volume  spectrum,  we  derive  a  2-D  horizontal  spectrum  through 
integration  yielding 


wjff  (K, L)  =  J  WjV) (K,L,M)dM  . 

(12) 

Substituting  Eq.(l  1)  into  Eq.(12),  gives  us 

W^VH\K,L)=^A  B\\\2(k2 +L2)  +  M2l\^~l  dM  . 

(13) 

n  o  J 

For  P  =  2  ,  Eq.  (13)  can  be  reduced  to 

<>(K,Z)=A^[A=(r+I2)p. 

(14) 

Also,  if  the  values  of  f?~5xl0  4  and  A~5  are  chosen  as  in  Yamamoto’s  work,'-12]  Eq.  (14) 

can  be  further  reduced  to 

wS{K,L)  =  a(K2+L2)-\ 

(15) 

where 

a  =  —  =  5xlCT5. 

(16) 

2A 

For  the  forward  propagation,  only  the  2-D  vertical  spectrum  in 

the  (x,  z)  plane  is 

required.  It  can  be  similarly  defined  by  reducing  the  spectrum  according  to 
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w!r) 


(K,M)=  J  w}y\K,L,M)dL  =  ^-^-  ^A2(K2  +  L2)+M 


2V 


2n 


(/?+ 2) 


dL.  (17) 


For  jB  =  2  ,  Eq.  (17)  reduces  to 


_  3 

wip(K,M)  =  a[25K2+M2Y 


(18) 


where  a'  =  1 .25  x  1 0  : . 

2.  Density  Fluctuations 

Many  of  the  inhomogeneities  within  the  bottom  volume  that  result  in 
reverberation  are  caused  by  variations  in  the  bottom  volume  density,  p.  The  effect  of 
these  variations  in  density  is  treated  though  the  definition  of  the  effective  index  of 
refraction, [8] 


'2  2  . 
n  =  n  +- 


2k: 


1  V72  3 

—Wp-~ 
P  2 


-Vpl 

P  ) 


(19) 


By  assuming  that  the  environment  is  approximately  range-independent  for  any  particular 
scattering  patch  and  that  sediment  properties  are  nearly  horizontally  stratified,  we  may 
simplify  this  expression  for  the  effective  index  of  refraction  to 


n  2  .  i 
n  =  n  -I - - 

2k: 


1  d2p 

p  dz2 


r  ]_dp  ^ 

p  dz 


(20) 


For  forward  propagation,  the  sound  speed  index  of  refraction  is  based  only  on 
large  scale  features,  such  that 


n;  ->  n2b(r,z ) 


(21) 


cb  is  defined  as 

c„=cbo{l  +  bz  +  Sl)  =  cbo  +  Sch  , 


(22) 
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where  cA  is  the  mean  bottom  sound  speed  at  the  interface,  b  is  the  normalized  gradient 

of  bottom  sound  speed  defined  by  bz  =  "/  ,  8 \  is  the  zero-mean  random  perturbation 

/  CK 

for  the  long  wavelength  component,  and  8ch  =  cK  (bz  +  Sl ) ,. 

Yamamoto’s  analysis  1  l2J  states  that  the  relative  fluctuations  in  density  relate  to  the 
relative  fluctuations  in  sound  speed  according  to  the  relationship 


Sp  _  2(pr-p0)  Sc  _  Sc 

Po  2  A)  “  Pr  C0  C0 


(23) 


where 


co  =c„0(l  +  bz)  and  y 


Pr-Po 

2PO~Pr 


(24) 


In  this  case,  p0  and  c0  are  the  ambient  values  of  the  density  and  sound  speed  in  the 

sediment,  respectively,  and  p  is  the  density  of  the  grain  (chosen  as  2650  3  ).  This 

/  m 

allows  us  to  express  the  bottom  volume  density  as 


P  Po  ^p  p0 


1  +  ^Sc 


(25) 


Taking  the  first  and  second  partial  derivatives  of  Eq.  (25)  with  respect  to  depth,  z,  and 
neglecting  depth  gradients  in  either  c0  or  p(] ,  we  obtain  expressions  for  both  the  change 
in  density  and  rate  of  change  in  density  with  depth 

=  (26) 

az  c0  oz 

and 


Ip_ 

dz2 


lypo  d2 

c0  dz2 


(Sc). 


(27) 


D.  TURBULENT  SOUND  SPEED  PERTURBATION 

Understanding  the  mechanism  of  reverberation  from  the  elements  of  the  bottom  is 
an  important  step  in  predicting  the  expected  reverberation  from  the  sea  floor.  However, 
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these  acoustic  signals  must  still  travel  through  the  water  column  on  their  way  from  source 
to  bottom  scattering  patch  and  from  the  scattering  patch  to  the  receiver.  Small-scale 
fluctuations  in  sound  speed  within  the  water  column  can  alter  the  paths  of  travel  for  these 
acoustic  signals  on  their  way  to  and  from  each  scattering  patch  and  therefore  affect  the 
overall  structure  of  the  reverberation.  Such  fluctuations  are  common  and  are  often  the 
result  of  oceanic  phenomena  such  as  turbulence  or  internal  waves. [17] 

Of  the  phenomena  that  generate  water  volume  sound  speed  fluctuations,  we  will 
place  our  emphasis  on  the  effect  of  turbulence  since  the  focus  of  our  later  analysis  will  be 
on  reverberation  signal  coherence.  To  characterize  the  turbulent  perturbation  field,  a 
statistical  model  is  used.  In  this  method,  the  turbulent  field  is  approximated  as  a  random 
realization  of  perturbations  based  on  an  accepted  power  spectral  density  of  the 
turbulence.  This  approach  should  be  consistent  with  perturbation  theory^,  which  is  only 
valid  for  the  weak  scattering  regime,  the  case  where  the  induced  change  in  the  signal 
strength  is  not  large  enough  to  exceed  the  average  signal  strength.[17] 

The  index  of  refraction  for  perturbations,  //(f) ,  is  defined  as 

Sn(r)  =  ^-  =  jU(r),  (28) 

c(r) 


where  Sc{r )  is  the  change  in  sound  speed  and  c(f )  is  the  sound  speed  along  the  path.  In 
this  case,  the  variance  of  //(f)  can  be  written  as 


var (//(f))  =  — f  S  (k)d3k , 

{271)  1° 


(29) 


where  S  (k)  is  the  spectral  density  of  the  perturbations  in  units  of  m  and  the  integral  is 
calculated  over  the  entirety  of  three  dimensional  (3-D)  wavenumber  space, 
k  =(K,L,M) .  Because  the  perturbation  region  is  not  bounded,  the  variance  is 


nonnalized  according  to 


var (//(f))  =  (|//(f)|2)  =  — - — f  /r(f)//V  (30) 

\  /  volume JV 

Combining  Eqs.  (29)  and  (30)  yields  the  relationship 


14 


\y^r^\y^k. 


(31) 


The  transfonn  of  the  index  of  refraction,  //(f),  to  the  signal  spectrum,  F  (k),  can  now  be 
defined  as 


Mir)  =  77-T  f  F(k)e,k*r  dkk  . 
\2n) 


(32) 


Using  this  transfonn,  the  correlation  of  ju(r )  can  be  written  as 


j  ju(r)ju\r  -r)d3r  =  j  d3r  j  FM(k)elk*r  dkk  j  F*t(k')e  lk*(r~r  ]dkk 

yZ7T)  _cx 


(33) 


where 


Fu{k) 


is  the  Fourier  transfonn  of  the  correlation  function.  The  intenelation 


between  the  square  of  the  signal  spectra,  F  (k),  and  the  spectral  density,  Sfl(k) ,  can  be 
determined  by  setting  r  to  zero  to  yield 


Fu  (*)  =  ( volume )  *  Su  (k) , 


which  can  be  rewritten  as 


F  (k)  =  volume* S  (k) 


n  1/2 


(34) 


(35) 


A  similar  process  holds  in  2-D.  Holding  one  of  the  variables  at  a  fixed  value  can 
create  a  planar  slice.  Holding  y  constant  at  y=0  yields 


- —  [[  l T(x,0,z)dxdz  = — [  [  V  (K,M)dKdM 

ron  (27l)  ^  ^  ^ 


Area 


(36) 


where  V  (K,M)  is  the  two-dimensional  (2-D)  spectral  density  along  the  plane  y  =  0. 
Furthermore,  we  can  also  represent  the  index  of  refraction  as 


=  f  GM(K,M)eHKx+Mz)dKdM , 

(2k) 


(37) 


15 


where  GU(K,M)  is  the  spectral  transform  of  the  2-D  perturbation  field.  The  fonns  of  Eq. 
(36)  and  (37)  suggest  a  relation  similar  to  (35), 

| Gm  ( K , M) |  =  [  Area  *  VM  ( K ,  M)]  1/2 .  (38) 


The  2-D  spectral  density,  VU(K,M),  can  also  be  derived  from  the  3-D  spectral  density, 
SM(k) ,  as [14] 


=  SM(K,L,M)dL.  (39) 

Now  that  we  have  developed  expressions  to  relate  sound  speed  fluctuations  to 
spectral  densities,  we  must  leverage  these  relationships  to  produce  the  turbulent 
perturbation  field.  This  field  can  be  approximated  through  the  use  of  the  method  of 
structure  functions. [  141  This  method  uses  both  isotropic  and  homogeneous 
approximations.  It  also  reduces  the  scale  of  the  analysis;  therefore  large-scale 
perturbations  are  not  accounted  for.  Using  standard  statistical  relations,  the  3-D  spectral 
density  function  can  be  written  in  the  form 1171 

E{k)  =  Aeink~xm ,  (40) 

where  fs(£)is  the  spectral  density,  A  is  a  scalar  multiplier,  £  is  the  energy  dissipation 

rate,  and  k 2  =  K2  +L2  +M 2 .  A  lower  wavenumber  threshold  limit,  kt.  is  also  established 
such  that  the  K 1/3  dependence  becomes 


E(k)~ 


1 

(k2+k;)W6  ' 


(41) 


This  prevents  E(k)  from  blowing  up  at  small  values  of  k.  A  highpass  filter  of  the  form 


k 2 

k2  +  kf 


(42) 


is  also  applied  to  force  E(k)=  0  when  k= 0.  A  high-frequency  cut-off  based  on  the 
Batchelor  spectrum  is  then  introduced  defined  by [17] 
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Rb  ( k )  =  exp 


-q 


\kB  j 


(43) 


where  q  is  an  order  unity  factor  (in  this  case,  a  value  of  3.7  will  be  used  to  maintain 
consistency  with  previous  work),[17] 

Batchelor  wavenumber  can,  in  turn,  be  expressed  as 


and  kB  is  the  Batchelor  wavenumber.  This 


(  ~  V/4 

£ 


kvkt  j 


(44) 


where  £  is  the  depth  averaged  kinetic  energy  dissipation  rate  in  units  of  W/kg,  v  is  the 
kinematic  viscosity  (1.40e-6  m /s),  and  kt  is  the  thennal  diffusivity  of  sea  water 
(~1.00e-7  m2/s). 

An  expression  for  the  depth  averaged  kinetic  energy  dissipation  rate  has  been 
developed  for  the  turbulence  caused  by  the  decay  of  internal  waves.  This  expression  was 
originally  developed  for  use  in  deep-water  problems  and  its  applicability  to  shallow  water 
remains  unknown.  In  this  development,  the  depth  averaged  kinetic  energy  dissipation 

rate,  £ ,  is  defined  as [15] 


f  =  0.3 


1 

7 


\ 4 

(  —j*N  (z)bEGMfj 

2 

cosh-1 

r  N(zj 7 

V  fl  )j 

(45) 


where  f  ]  is  the  inertial  frequency  in  units  of  rad/sec,  y'»  is  the  mode  number  (  j,  =1  is 
the  characteristic  mode  for  shallow  water),  N(z)  is  the  buoyancy  frequency  in  rad/sec,  and 
bEGM  is  a  measure  of  internal  wave  intensity  (=0.5 m).  The  notation  (  ^  denotes  a 
depth  average.  The  buoyancy  frequency,  N(z),  can  be  computed  using 


N(z)  =  gin 


1 

P(z) 


(46) 


where  p(z)  and  p  (z)  are  the  density  and  the  potential  density  as  a  functions  of  depth, 
respectively. 
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Following  this  development,  the  3-D  spectral  density  can  now  be  written  in  the 


form 


SM(k;z) 


AAz) 

(k2+k2)W6 


(47) 


where  ^4  (z)is  introduced  as  the  turbulent  strength  parameter  in  units  of  m' /3J1J|  This 
parameter  is  defined  as 


AM(Z)  =  2.6a2rm 


—j,N(z)bEGM 

\7T 


y/3  , 

V  dT'^ 

2\ 

( ln(2iV(z)//7)^ 

2  \ 

^  dz  J 

/ 

(2 N(z)/f,)  ) 

(48) 


where  a  =  S^e-S^C)'1,  rm=0.2,  and 


'ar'V 

V  J 


is  the  variance  of  the  temperature 


gradient.  This  parameter  varies  slowly  with  respect  to  the  correlation  length  scales  of  the 
turbulence,  so  the  previous  isotropic  assumption  holds  locally. 

To  develop  a  2-D  spectral  density,  integration  can  be  applied  to  the  portion  that 
varies  with  the  turbulent  spectrum  only.  The  remaining  terms  such  as  the  turbulent 
strength  parameter  and  wavenumber  roll-off  can  then  be  included  afterwards.  The 
portion  of  spectral  density  that  varies  with  the  turbulent  spectrum  only  can  be  written  as 


*3  n(k) 


1 


{k-+k;) 


2 \1 1/6 


The  2-D  spectral  density  is  derived  through  integration  as 


(49) 


<$>2D{K,M)  =  -1-  f  d >3D(K,L,M)dL 
in  J 


-J 

J  TT  J 


1 


2n  (K  +  L  +  M  +  k') 


2x11/6 


<IL . 


(50) 


Substituting  a  =  K2  +M1  +  k2  into  Eq.  (50)  yields 
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<t \D(K,M) 


— !_f 

nau /3 { 


(1+i/ 


a’ 


11/6 


dL . 


(51) 


Solving  this  integral  (from  an  integral  table)  produces  the  expression 

i  V^r(|) 


®2D(K,M). 


2nam  IW) 


(52) 


When  the  terms  that  do  not  vary  with  the  turbulent  spectrum  are  reapplied,  the  resulting 
expression  for  2-D  spectral  density  of  turbulent  fluctuations  becomes 


VM(k',z) 


V^r(j)  am(z)  e 

In  r(^)  (k2+kf)4'3  k2  +  kf 


(53) 


E.  CORRELATION  METHOD  OF  SEPARATING  AND  DISCRIMINATING 

VOLUME  AND  SURFACE  SCATTERING 

In  the  shallow  water  environment,  the  primary  contributors  to  reverberation  are 
the  roughness  of  the  bottom  interfaces  and  the  inhomogenieties  within  the  bottom 
volume.  A  reverberant  return  in  this  type  of  environment  is  essentially  the  combined 
effect  of  these  two  factors.  In  some  cases,  it  would  be  desirable  to  have  a  metric  for  the 
degree  to  which  each  factor  contributes  to  the  total  reverberation.  For  example,  methods 
for  determining  the  bottom  sediment  type  and  extracting  certain  bottom  parameters  from 
the  reverberant  return  have  been  proposed. [13]  However,  these  methods  require  that  the 
primary  contributor  to  the  reverberation  be  known.  The  influence  of  changes  in  source 
and  environmental  parameters  on  the  reverberation  contributions  from  both  the  interfaces 
and  the  bottom  would  be  useful  as  well. 

To  separate  the  contributions  of  the  interface  roughness  and  volume  perturbation 
components,  the  following  approach  is  suggested  based  on  the  work  of  Ivakin.[13]  This 
approach  involves  the  measurement  of  the  spatial  coherence  of  the  reverberant  field.  If 
the  two  contributors  are  uncorrelated  with  respect  to  each  other,  then  the  contribution 
from  interface  roughness,  pr,  and  the  contribution  from  bottom  volume  perturbations,  pv, 
have  intensities,  /,-  and  Iv  respectively,  that  are  additive  as 
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=Ir+Iv, 


(54) 


1  = 


where  /  is  the  total  intensity  and  p=  pr+  pv,  the  coherent  sum  of  the  two  reverberation 
pressure  signals.  The  spatial  coherence  coefficient,  N,  can  then  be  defined  as 


N(rx,ri) 


(p{rx)p’{r2)) 


(55) 


where  rj  and  r2  are  the  spatial  positions  of  two  receiving  elements.  Once  the  distance 
between  the  receiving  elements  and  the  scattering  patch,  R ,  becomes  sufficient  such  that 
L  =  /,  +  r\<$cR,  the  resulting  spatial  coherence  coefficient  becomes 


N(L)=j-N,.(L)  +  j-N,(L). 


(56) 


If  these  partial  roughness  and  volume  coherence  coefficients,  Nr  and  Nv  for  the  interface 
roughness  and  bottom  volume  perturbations,  respectively,  have  significantly  different 
coherence  scales,  then  the  contribution  of  each  component  can  be  detennined  through 
measurements  of  the  spatial  coherence  of  the  reverberant  field. 

Consider  the  case  of  a  near  backscattering  geometry  in  which  the  bottom  volume 
has  the  same  sound  speed  as  the  water  and  the  attenuation  coefficient  (3  =  2k  S  where  S  is 
the  loss  parameter  for  a  particular  bottom  type.  The  expression  for  the  partial  roughness 
coherence  coefficient  can  be  obtained [13]  as 


Nr(L)  =  eikL^x 


in  (L/Lr) 

{L/Lr)  ’ 


(57) 


2H(T/t) 

where  X  is  the  grazing  angle  and  Lr  = - ? - for  a  period,  T,  and  a  duration,  T  .  This 

;rsin  / 

approximation  holds  for  the  case  when  T 12  «:  R/c .  A  similar  expression  is  available  for 
the  partial  volume  coherence  coefficient  in  the  form  of 

1 


(V,(£)  = 


_  ikLsinx 


1  ~i{L/Lv) 


(58) 
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where  Lv  = 


4H5 
sin2  x 


This  approximation  holds  for  the  case  when  1  /  jB  <sc  R  . 


For  this  case,  the  expression  for  the  partial  roughness  coherence  coefficient,  Nr, 
displays  a  dependence  on  the  signal  duration,  x ,  while  the  expression  for  partial  volume 
coherence  coefficient,  Nv,  displays  a  dependence  on  the  sediment  parameters  by  way  of 
the  loss  parameter,  S .  The  different  effects  of  the  two  partial  coherence  coefficients  can 
be  more  easily  seen  if  parameters  are  chosen.  The  cases  where  H=30m,  ^=60°, 
=0.005  (which  corresponds  to  the  loss  parameter  for  silt),  and  signal  durations,  x ,  of 
107  and  5 T  developed  in  Ivakin’s  work  '■13]  are  shown  in  Figure  1.  In  this  figure,  the 
different  natures  of  the  partial  roughness  coherence  coefficients  are  clearly  visible. 


(a)  (b) 

Figure  1.  Coherence  Coefficient  Magnitude  Plots  for  Two  Pulse  Durations,  (a)  x=\0T 

and  (b)  x  =5T  (from  [13]) 


This  method  describes  an  interesting  and  potentially  powerful  analysis  tool.  It 
also  suggests  a  similar  frequency-domain  analog  for  broadband  signals.  In  the  previous 
example,  pulse  lengths  are  chosen  with  respect  to  period  such  that  x IT  is  held  constant 
at  values  of  10  and  5,  respectively.  This  is  analogous  to  holding  the  ratio  between  center 
frequency,^,  and  bandwidth,  B W,  fixed  at  the  same  constant  such  that  x IT  =  fcl BW  . 
As  such,  similar  differences  in  spatial  coherence  should  be  produced  by  employing 
different  bandwidths,  produced  either  by  varying  the  source  spectrum  or  band-pass 
filtering  of  the  received  broadband  signal.  Given  an  accurate  model  to  predict  the 
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expected  spatial  coherence  for  a  given  environment  and  geometry,  spatial  coherence 
analysis  of  broadband  data  analyzed  at  different  bandwidths  produces  consistent  results 
as  the  narrow-band,  pulse  length  varied  theory  indicates. [19] 
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III.  NUMERICAL  METHODS 


A.  MONTERE Y-MIAMI  PARABOLIC  EQUATION  (MMPE)  MODEL 

As  the  name  suggests,  the  MMPE  Model181  is  based  on  the  parabolic 
approximation  to  the  acoustic  wave  equation.  This  approach  begins  by  first  representing 
the  time  hannonic  acoustic  field  in  cylindrical  coordinates.  Azimuthal  symmetry  is  then 
assumed  since  ocean  environments  tend  to  display  weak  azimuthal  dependence.  The 
resulting  acoustic  field  can  then  be  represented  as 

P{r,z,cot )  =  p[r,z)e~,m  .  (59) 


Substituting  Eq.  (59)  into  the  wave  equation  in  cylindrical  coordinates  yields  the 
Helmholtz  equation, 

V2p(r,z)+  ^  p(r ,z)  =  0 ,  (60) 

c(r,zf 

where 


i  a  a  a2 

- r - 1 - - . 

r  dr  dr  dz 


(61) 


The  Helmholtz  equation  can  then  be  factored  by  introducing  an  operator,  Qop ,  with  the 
notation 

Qop=(ji  +  e  +  l)K  (62) 

where 

i  a2  2  , 

A  =  7TTT,  £  =  n  “I.  (63) 

k-  dz1 


Q 

and  the  index  of  refraction  n  =  — ,  where  Co  is  the  reference  sound  speed  of  the  ocean 

c 

volume. 

Once  the  effect  of  cylindrical  spreading  is  considered  and  the  Helmholtz  equation 
is  properly  factored,  the  out-going  acoustic  pressure  can  then  be  defined  as [9], 
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(64) 


P{r’z)  =  Po^Q~p2w{r,z)e‘Kr , 

where  y/{r,z)  is  the  envelope  function  or  PE  field  function  and  P0  is  the  pressure 
amplitude  at  Ro .  The  parabolic  equation  for  the  PE  field  function,  yAj~,z),  is  then  defined 
by 

^7  =  ~iko¥  +  iKQopV  =  ~ikoHoP¥ ,  (65) 

where  Hop  is  a  Hamiltonian-like  operator  defined  by 

Hop=l~Qop,  (66) 

which  defines  the  evolution  of  the  PE  field  function  in  range.  The  relationship  between 
values  of  ^at  different  ranges  can  be  defined  by 

l//(r  +  Ar)  =  <J>(r)^(r) ,  (67) 


where  <J>(r)  is  a  propagator  that  accounts  for  the  change  in  the  values  of  y/ with  range.  A 
split-step  Fourier  (PE/SSF)  method1101  can  then  be  used  to  determine  this  propagator, 
<f>(r).  The  Fast  Fourier  Transform  (FFT)  allows  the  PE/SSF  implementation  to  be 
represented  by[8] 


l//(r  +  Ar,z)  =  e 


■I  Al'r, 

-,k"Yu« 


(r+Ar.z) 


FFT 


-ik0ArT  (k. ) 


IFFT 


-iK^op(r,z) 

e  W\r,z) 


V 


(68) 


where 


and 


T  =1- 

op 


1- 


rk_  ^ 

V  ko  J 


-il/2 


(69) 


(70) 
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B.  NUMERICAL  IMPLEMENTATION 

At  this  point,  we’ve  described  the  applicable  concepts  of  reverberation  theory, 
particularly  the  effects  of  scatter  from  bottom  interfaces  and  the  bottom  volume  in  terms 
of  reverberation  loss.  We’ve  also  described  methods  for  generating  perturbations  in  the 
interface  roughness,  the  bottom  volume,  and  the  turbulent  sound  speed  in  the  water 
column.  The  next  step  is  to  apply  all  of  this  theory  in  such  a  manner  as  to  produce  useful 
results  from  numerical  analysis.  The  following  describes  the  application  of  this  theory  to 
the  MMPE  Model. 

1.  Interface  Roughness 

In  the  discussion  of  interface  perturbation  theory,  the  1-D  interface  spectrum  was 
developed  (Eq.(9)).  In  order  to  generate  a  1-D  roughness  realization  from  this  1-D 
spectrum,  we  define  the  roughness  realization  as 

7](x)  =  j^St(K)eXxdK,  (71) 

where 

S, ( K )  =  [w"  ( K)Y2  A{K)emK) .  (72) 

A(K)  provides  a  random  amplitude  and  9(K )  provides  a  random  phase.  Combined, 

these  components  form  a  randomization  factor  of  the  form,  A  e'e ,  for  which  the  complex 
amplitude  exhibits  a  normal  distribution.  To  accomplish  this,  the  phase  component  is  a 
random  variable  that  is  distributed  uniformly  over  the  interval  [0,2;r].  In  practice,  this 
component  is  generated  by 

9  =  2m\  ,  (73) 

where  rx  is  a  uniformly  distributed  random  variable  over  the  interval  [0,1].  The 
amplitude  component  is  a  random  variable  with  a  Rayleigh  distribution  generated  by 


where  r2  is  another  unifonnly  distributed  random  variable  over  the  interval  [0,1].  When 
applied  to  the  MMPE  model,  the  actual  relationship  used  is 
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(75) 


<to=(i+4,KT” 

nix) 

and  the  final  result  is  normalized,  i.e.  — — — - .  The  actual  rms  perturbation  is  then 

(n2(x))2 

achieved  by  simply  multiplying  the  result  by  the  desired  rms  value. 

2.  Bottom  Volume  Sound  Speed  Perturbations 

In  the  discussion  of  bottom  volume  perturbation  theory,  the  2-D  vertical  spectrum 
was  developed  (Eq.  (17)  and  (18)).  In  a  fashion  similar  to  that  used  for  1-D  interface 
roughness  realizations,  a  2-D  vertical  volume  sound  speed  fluctuation  realization  is 
defined  as 


c0S(x,z)  =  J  J  S2{K,M)eiKxeiMzdKdM  , 


(76) 


where 


S2(K,M ): 


W2P  (K,M) 


A(K,M)e 


iO{K,M ) 


(77) 


Both  A(K,M)  and  6{K,M )  are  similar  to  the  random  values  A(K )  and  0{K)  introduced 
in  Eq.  (72).  Specifically,  they  are  generated  by 

0{K,M)  =  2nrx{K,M)  (78) 

and 

A(K,M)  =  ,/-ln (r2(K,M))  ,  (79) 


respectively,  where  both  )  and  r2{K,M)  are  now  a  matrix  of  unifonnly 

distributed  random  numbers  in  [0,1].  When  applied  to  the  MMPE  model,  the  actual 
relationship  used  is 

W^\k,M)oc(\2K2 ,  (80) 

and  the  result  is  normalized  and  rescaled  by  the  desired  rms  values. 
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3.  Bottom  Volume  Density  Fluctuations 

The  definition  for  2-D  vertical  volume  sound  speed  fluctuation  realizations  (Eq. 
(76))  can  also  be  combined  with  the  results  from  density  fluctuation  theory  (Eqs.  (26)  and 
(27))  to  provide  a  realization  for  the  bottom  density  fluctuations  as  well.  Specifically, 
substituting  Eq.  (76)  into  Eqs.  (26)  and  (27)  yields  both 

[  [  (iM)S2  (K,M) e,Kx eiMz clKclM  (8 1 ) 

az  c0 

and 

^£_  =  ^]PiLl  f  (-m2)S2  (K,M)e,KxelMz  clKclM  .  (82) 

az  c0 


The  sediment  effective  index  of  refraction  is  then  derived  by  substituting  Eq.  (81)  and  Eq. 
(82)  into  the  expression  for  the  effective  index  of  refraction,  n'2 ,  Eq.  (20),  yielding 

n2  (x,z)  =  n l  (x,z)  +  -^y[a(x,z)  +  /3{x,z)\  ,  (83) 

2k0 

where 


a(x,z)=  2r  ff M2S,  (K,M)e/KxeiM:dKdM 
c0(x,z)JJ 


(84) 


and 


jB(x,z)  =  + 


=  4-3/ 


2/ 


-12 


,{x,z) 


JJ MS2  (K,M  )e,KxeMzdKdM 


(85) 


Both  of  these  parameters  are  computed  in  parallel  with  8ch  during  the  implementation  of 
the  MMPE  model.  An  additional  propagator  term  is  then  used  in  the  sediment 

<t>p{x,z)  =  eiArk°uo{x’z),  (86) 

where 

^p(x’z)  =  _7r7y[a(x’z)  +  ^(x’z)]-  (87) 

4/c0 
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4.  Water  Volume  Sound  Speed  Turbulence 

The  relationship  between  the  2-D  spectral  density  of  turbulent  fluctuations, 
V  (k :;z)  ,  and  the  2-D  index  of  refraction  for  perturbations,  ju(x,z)  ,  (also  known  as  the 

perturbation  field)  is  evident  in  Eqs.  (37)  and  (38).  Using  this  relationship,  random 
realizations  of  ju(x,  z )  are  generated  by  modifying  Eq.  (37)  to  create  the  expression 

ju{x,z)  =  j  j  eKKx+Mz)dKdM ,  (88) 

where  A(K,M)  provides  a  2-D  random  amplitude  for  each  wavenumber  component  and 
0(K,M )  provides  a  2-D  random  phase  component  similar  to  the  random  values 
A(K,M)  and  0(K,M )  introduced  in  Eqs.  (79)  and  (80). 

The  2-D  sound  speed  perturbation  field  is  then  easily  obtained  from  this 
realization  of  the  perturbation  field  by  multiplying  by  ca ,  since  8c(x,  z)  =  c0ju(x,  z) .  At 

this  point,  randomizations  in  the  perturbation  field  have  been  generated,  but  no  regard  has 
been  given  to  scale.  As  before,  after  the  realization  of  the  perturbation  field  is  created, 
the  result  is  simply  nonnalized  and  rescaled  by  the  desired  rms  value. 

C.  TIME-DOMAIN  PROCESSING 

The  preceding  development  has  covered  reverberation  loss  theory  and  developed 
a  method  for  its  numerical  implementation.  This  development  is  based  on  continuous 
wave  (CW)  analysis,  which  holds  for  discrete  frequencies.  However,  to  obtain  the 
structure  of  the  multipath  pulse  propagation  in  time,  broadband  analysis  is  required.  To 
create  a  broadband  source,  the  MMPE  model  is  run  for  a  large  number  of  discrete 
frequencies  to  generate  an  output  over  a  given  frequency  band.  Time  domain  analysis 
then  extracts  the  nature  of  the  two-way  travel  time  structure  of  the  reverberation  loss 
from  both  the  interface  and  the  volume  reverberation.  The  reverberant  field  at  each  range 
step  can  then  be  detennined  and  the  propagation  continued  through  the  water  column. 

1.  Time-Domain  Analysis  of  Interface 

The  basic  geometry  of  a  typical  two-way  return  from  a  scattering  patch  as  adapted 
from  Smith  and  Cushman1^1  is  depicted  in  Figure  2.  The  travel  time  structure  of  the  one¬ 
way  forward  propagating  field  at  a  given  range  step  rm  is  denoted  as  p+(rm,z,t),  where  t  is 
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r 

the  reduced  time  defined  as  t  =  T  -  —  ,  and  T  is  the  actual  travel  time.  The  two-way 

co 

pressure  field  as  observed  at  the  receiving  element  is  determined  by  performing  the 
convolution  of  two,  one-way  propagated  fields  in  the  time-domain^41 

P2-Way,b  (r,„  >0  =  {  P+Tb  (>;,  >  0 P+Rb  (rm  ’  T  ~  t)dt  ,  (89) 

where  p+rb (rm , t)  and  plRb(rm,T-t)  are  the  forward  propagated  pressure  fields  from  the 
transmitter  and  receiver  to  the  scattering  point  as  evaluated  at  the  bottom  interface  for  the 
range  step  rm,  respectively,  and  ris  the  reduced  time  of  the  reverberation  defined  as  T=  T 
It 

- — .  Reciprocity  allows  the  forward  propagating  pressure  field  from  receiver  to 

co 

scattering  patch,  p+Rb ,  to  replace  the  return  path,  p_Rb ,  since  these  two  propagated  fields 
are  identical. 


Multi-Static  Reverberation  Geometry 


Figure  2.  Bistatic  Return  From  Scattering  Patch 
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In  practice,  this  time-domain  convolution  of  the  two  pressure  fields  can  be  more 
easily  calculated  in  the  frequency  domain,  where  the  operator  becomes  a  scalar 
multiplication  of  the  two  functions.  This  allows  the  two-way  field  from  the  interface  to 
be  expressed  in  the  frequency  domain  as 

Pl-.ay.b  (f„  >  /)  =  P +Tb  fa, ,  f)P+Rb  fa ,  /)  ,  (90) 

where 

P+Tb  fa ,  /)  =  -T=  V+Tb  ( rm >  f)^  (9  1 ) 

Vr». 

and 

P+Rh  fa > /)  =  -7=  V+Rh  fa > />'V  •  (92) 

"V  f » 

For  each  bottom  scattering  patch,  the  two-way  travel  time  structure  of  the  reverberation 
loss  for  the  bottom  interface,  RLb ,  can  then  be  defined  as 

P-b(rm’t)  =  A\ P2-way,b(rn,J)e~i27r,'df  ’  (93) 

where  the  constant  A  is  included  to  incorporate  all  the  other  factors  needed  to  define 
reverberation  loss,  RLb.  This  provides  the  two-way  travel  time  structure  due  to  scattering 
from  a  single  range  step,  rm.  This  calculation  is  then  conducted  for  each  range  step,  and 
the  total  field  at  the  receiver  is  computed  as  the  coherent  sum  of  the  pressure  values  from 
all  of  the  range  steps,  rm.  This  is  accomplished  by  matching  up  the  discrete  arrival  times, 
tn,  according  to  the  relation 

M 

P-b{tn)  =  YjP-b(rmAn)i  (94) 

m= 1 

where  p.}}  is  the  total  interface  reverberation  pressure  received  at  the  receiver  at  time  tn. 

2.  Time-Domain  Analysis  of  Volume 

The  time  domain  analysis  is  more  complicated  for  the  bottom  volume  since 
multiple  depth  steps  must  also  be  taken  into  account.  In  a  manner  similar  to  that  used  for 
the  interface  analysis  (Eq.  (90)),  the  reverberant  field  due  to  each  point  in  range  and 
depth  is  computed  by  combining  the  source-to-patch  and  the  patch-to-receiver 

propagating  field  in  the  frequency  domain  using 
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P  2-way  ,v  ( rm  >ZJ)  =  nK  ,  z)p+T  fa  ,  Z,  f)p+R  fa,  ,Z,f), 


(95) 


where  rm  is  the  range  step,  z  is  the  depth  step  ( z  >  zh  for  the  bottom  volume),  /  is  the 
frequency,  and  nfa,z)  is  the  local  index  of  refraction.  A  Fourier  transform  yields  the 
appropriate  time-domain  response 

P2-Way,Vfa’Zfa  =  B\ >  (96) 

where  B  is  a  constant  that  accounts  for  all  the  other  terms  needed  to  define  the 
reverberation  loss,  RLV.  The  two-way  travel  time  structure  of  the  reverberation  loss  for 
the  volume,  RLV,  can  then  be  derived  from 

P-vfa>t)=  J  P2-Way,vfa’Z’t)dz  ’  (97) 

z>zb 

which  is  the  coherent  sum  over  all  depths  below  the  interface  at  a  range  step  m.  Each 
single  set  of  time  series  can  then  be  matched  and  summed  as 

M 

P-V(t«)  =  'ZP-v(rm’tn)’  (9g) 

m= 1 

where  p.v  is  the  two-way  time  domain  pressure  defining  the  volume  reverberation  loss  as 
measured  at  the  receiver  at  time  t„  due  to  the  entire  bottom  volume. 

D.  MODELING  PARAMETERS 

1.  Numerical  Analysis  Geometry 

The  use  of  a  vertical  linear  array  (VLA)  during  ASIAEX  allowed  for  the  capture 
of  varying  spatial  properties  in  the  acoustic  return  along  its  axis.  The  utility  of  analyzing 
such  spatial  variation  has  been  shown  through  a  discussion  of  the  correlation  method  of 
separating  and  discriminating  volume  and  surface  scattering.1 [13]  However,  this  method 
requires  an  accurate  prediction  of  the  expected  spatial  response  of  the  interfaces  and  the 
bottom  volume  for  a  given  geometry.  It  is  therefore  desirable  to  generate  a  similar 
vertical  array  for  the  numerical  analysis  to  produce  the  required  prediction.  In  this  case,  a 
16-element  array  was  modeled.  This  modeled  VLA  extends  60m,  from  20m  to  80m  in 
depth,  with  a  constant  4m  vertical  separation  between  elements  as  shown  in  Figure  3. 
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Multi-Static  Reverberation  Geometry 


Figure  3.  Geometry  of  VLA  and  Scattering  Patch 

The  problem  geometry  and  environmental  parameters  are  chosen  to  closely 
simulate  the  conditions  observed  during  the  ASIAEX.  This  modeled  environment 
includes  both  a  water/bottom  interface  at  an  average  depth  of  100m  and  a  bottom/sub- 
bottom  interface  at  an  average  depth  of  105m  -  to  simulate  both  the  smoother,  softer  and 
more  penetrable  bottom  layer  of  sand  and  mud  -  and  a  harder  sub-bottom  layer  of  coarser 
sand,  gravel  or  rocks.  The  element  at  48m  is  chosen  as  the  source,  since  it  is  close  to  the 
set  charge  depth  of  50m.  The  sound  speed  profde  (SSP)  was  taken  directly  from  CTD 
measurements  taken  during  the  ASIAEX,  averaged  over  the  duration  of  the  exercise.  The 
choice  of  frequency  bands  and  the  geometry  of  the  vertical  array  geometry  are  made  to 
approximate  ASIAEX  as  close  as  possible  while  still  allowing  for  reasonable  calculation 
speeds.  The  maximum  propagation  range  is  5km  to  capture  most  of  the  environmental 
contributors  to  the  reverberation  loss. 
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To  ensure  that  the  results  of  the  numerical  analysis  hold  for  a  variety  of 
environmental  parameters,  the  numerical  analysis  was  also  conducted  for  a  number  of 
environmental  variations.  In  each  case,  only  one  parameter  was  varied  from  the 
“reference.”  The  parameters  for  this  reference  environment  are  a  lm  rms  variation  in  the 
water/bottom  interface  roughness,  a  2m  nns  variation  in  the  bottom/sub-bottom  interface 
roughness,  a  15  m/s  rms  bottom  volume  sound  speed  perturbation,  and  no  water  volume 
sound  speed  turbulence.  All  comparisons  of  environmental  parameters  are  done  with 
respect  to  a  single  realization  of  this  reference  model. 

2.  The  Modeled  Environment  and  Variations 

a.  Variations  in  Interface  Roughness 

Three  variations  in  bottom  interface  roughness  were  generated.  In 
addition  to  the  reference  condition  of  lm  water/bottom  interface  roughness  and  2m  rms 
bottom/sub-bottom  interface  roughness,  a  lower  roughness  realization  of  0.5m  and  lm 
nns  and  a  higher  roughness  realization  of  2m  and  4m  rms  water/bottom  and  bottom/sub- 
bottom  interface  roughness,  respectively,  were  calculated  and  analyzed.  Graphical 
representations  of  all  three  realizations  are  depicted  in  Figure  4. 

b.  Bottom  Volume  Sound  Speed  Perturbations 

To  model  an  environment  with  a  soft  bottom  and  a  harder  sub-bottom,  the 
sound  speed  of  the  bottom  layer  was  set  at  1700m/s  while  the  sound  speed  of  the  sub¬ 
bottom  layer  was  set  at  1760m/s.  Three  variations  in  bottom  volume  sound  speed 
perturbations  were  also  generated.  In  addition  to  the  reference  condition  of  a  15m/s  rms 
bottom  volume  sound  speed  perturbation,  a  lesser  perturbation  realization  of  5m/s  rms 
and  a  greater  perturbation  realization  of  45m/s  nns  bottom  volume  sound  speed 
perturbation  were  calculated  and  analyzed.  This  greater  perturbation  realization  is 
somewhat  more  dramatic  than  what  is  expected  to  occur  naturally,  but  it  is  still  useful  as 
a  means  of  comparison.  Graphical  representations  of  all  three  realizations  are  depicted  in 
Figure  5. 
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Figure  4.  Plots  of  water/bottom  and  bottom/sub-bottom  interface  profiles  for  increasing 
rms  roughness,  (top  right)  0.5m  &  lm;  (top  left)  lm  &  2m;  (bottom)  2m  &  4m  (water/ 
bottom  and  bottom/sub-bottom  rms  roughness  respectively) 


Figure  5.  Plots  of  environmental  sound  speed  for  increasing  bottom  volume  sound 
speed  perturbation,  (top  left)  5  m/s;  (top  right)  15  m/s;  (bottom)  45  m/s  rms  perturbation 
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c.  Water  Volume  Turbulence 

Five  variations  in  water  volume  sound  speed  turbulence  were  generated  as 
well.  In  this  case,  the  reference  condition  is  an  environment  with  no  water  volume  sound 
speed  turbulence.  The  four  additional  realizations  were  for  increasing  degrees  of 
turbulence,  specifically  0.5m/s,  lm/s,  2.5m/s,  and  5  m/s  rms  water  volume  sound  speed 
turbulence.  Graphical  representations  of  the  four  non-reference  realizations  are  depicted 
in  Figure  6. 
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Figure  6.  Plots  of  environmental  sound  speed  for  increasing  water  volume  sound  speed 
turbulence,  (top  left)  0.5  m/s;  (top  right)  1  m/s;  (bottom  left)  2.5  m/s;  (bottom  right)  5 

m/s  rms  turbulence 


3.  Source  Spectra 

The  output  of  the  MMPE  model  is  the  frequency  domain  response  of  the  field 

over  the  specified  bandwidth.  The  source  strength  for  each  frequency  bin  is  the  same, 

which  may  be  considered  like  a  square  spectrum.  Despite  the  fact  that  the  explosive 

sources  used  in  ASIAEX  are  wide  band  sources,  their  expected  frequency  spectrum  is  not 

flat.  The  actual  frequency  spectra  for  the  explosive  charges  used  in  ASIAEX  are  not 
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available,  but  the  explosive  spectrum  for  a  Standard  Underwater  Source  (SUS)  charge  for 
the  same  size  (38g)  and  detonation  depth  may  be  modeled  as  depicted  in  Figure  7.  The 
spectrum  for  a  1kg  SUS  charge  is  very  similar. 


Norm«iaed  SUS  Prettire  Spectnrn  wind o* 


Figure  7.  Normalized  Expected  SUS  Pressure  Amplitude  Spectrum  (at  lm  from  source) 

for  38g  charge  exploded  at  a  depth  of  48m 

In  previous  theses, [5]’ [6]’ [7]  a  Hanning  window  spectrum  was  applied  to  the  output 
of  the  MMPE  Model.  Though  this  spectrum  is  advantageous  from  a  signal  processing 
standpoint,  it  does  not  approximate  the  spectrum  expected  from  an  explosive  charge  well. 
Instead,  the  SUS  spectrum  over  the  calculated  bandwidth  is  applied  over  the  MMPE 
model  output.  This  is  conducted  as  a  multiplication  in  the  frequency  domain.  However, 
a  direct  comparison  with  previous  results  is  useful,  so  a  separate  set  of  calculations 
applies  a  Hanning  window  to  the  MMPE  Model  output  as  well. 

To  capture  a  large  portion  of  the  acoustic  energy  a  500  Hz  bandwidth  with  a  350 
Hz  center  frequency  calculated  over  1024  frequencies  was  chosen  for  the  analysis.  To 
approximate  an  explosive  source  as  used  in  ASIAEX,  the  spectrum  for  a  38g  charge 
exploded  at  48m  was  applied.  For  comparison  with  the  frequency  band  used  in  previous 
work,  a  250  Hz  bandwidth  with  a  250  Hz  center  frequency  calculated  over  512 
frequencies  was  also  analyzed.  To  make  a  closer  analysis  of  the  effect  of  changing  the 
center  frequency,  a  250  Hz  bandwidth  with  a  350  Hz  center  frequency  also  calculated 
over  512  frequencies  was  analyzed  as  well.  Both  the  SUS  charge  spectrum  and  a 
Hanning  window  were  applied  to  each  of  these  three  bands  for  a  total  of  six  variations. 
These  six  applied  spectra  are  depicted  in  Figures  8,  9,  and  10.  Only  the  500  Hz 
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bandwidth  SUS  spectrum  band  was  used  in  calculations  for  conditions  other  than  the 
reference. 


Frequency  iH*>  Frequency  iHz) 

Figure  8.  Hanning  and  Normalized  SUS  Spectrum  Windows  for  250  Hz  Bandwidth, 

250  Hz  Center  Frequency 


Frequency  (Hz*  Frequency  I  Mz) 

Figure  9.  Hanning  and  Normalized  SUS  Spectrum  Windows  for  250  Hz  Bandwidth, 

350  Hz  Center  Frequency 
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Figure  10.  Hanning  and  Normalized  SUS  Spectrum  Windows  for  500  Hz  Bandwidth,  350  Hz 

Center  Frequency 
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4.  MMPE  Input  Parameters 

Seven  different  input  files  are  used  to  define  the  desired  problem  parameters  when 
running  the  MMPE  Model.  The  infonnation  in  these  input  files  define  the  characteristics 
of  the  environment,  source,  receiver,  and  the  resolution  of  the  calculations.  The 
parameters  used  in  these  input  files  for  the  numerical  analysis  are  provided  in  Table  1 


Filename/Parameter 

Value 

Remarks 

Main  Control  File:  pefiles.inp 

Number  of  depth  points 

256 

Radix-2  integer  required  for  FFT 

Minimum  depth 

0  m 

Maximum  depth 

400  m 

Number  of  range  steps 

833 

Minimum  range 

0  m 

Maximum  range 

5.0  km 

Range  step  size 

6  m 

Maximum  computed  depth 

400  m 

Reference  sound  speed 

1500  m/s 

Source  File:  pesrc.inp 

Source  depths 

Varying 

Array  elements  at  20,  24,  28,  32,  36, 

40,  44,  48,  52,  56,  60,  64,  68,  72, 

76  and  80  m  depths. 

Center  frequency 

2  sets 

250  Hz  and  350  Hz 

Frequency  bandwidth 

2  sets 

250  Hz  used  for  both  input  center 
frequencies.  500  Hz  also  used  for 
350  Hz  center  frequency  in  a 
separate  run. 

No.  of  Frequencies 

2  sets 

1024  for  500  Hz  bandwidth,  512  for 
250  Hz  bandwidth.  Radix-2 
integer  required  for  FFT 

Sound  Speed  File:  pessp.inp 

Water  column  sound  speed 

2  sets 

SSPAvg  &  SSPPrevious 

Range  independent 

No.  of  SSPs  points 

2  sets 

84  for  SSPAvg  &  58  for  SSPPrevious 

Bathymetry:  pebath.inp 

Mean  bottom  depth 

100  m 

Range- independent  mean 

No.  of  depth  points 

1 

Bottom  properties:  pebotprop.inp 

Bottom  sound  speed 

1700  m/s 

Mean  value 

Sound  speed  gradient 

1  /s 

Mean  value 

Relative  density 

1.6 

No  density  variation 

Compressional  attenuation 

0.15  dB/km/Hz 

Shear  speed 

0 

Negligible 

Shear  attenuation 

0 

Negligible 
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Filename/Parameter 

Value 

Remarks 

Sub-Bottom  Bathymetry:  pedbath.inp 

Mean  Bottom  Depth 

105  m 

Range- independent  mean 

Sub-Bottom  Properties:  pedbotprop.inp 

Sub-bottom  sound  speed 

1760  m/s 

Mean  value 

Sound  speed  gradient 

1 

Mean  value 

Relative  density 

2 

Mean  value 

Compressional  attenuation 

0.2  dB/km/Hz 

Shear  speed 

150 

Constant 

Shear  attenuation 

0.5 

Constant 

RMS  Perturbations  (input  during  MMPE  run) 

Water/bottom  interface  roughness 

3  sets 

0.5m,  lm,  and  2m 

Bottom/ sub-bottom  interface  roughness 

3  sets 

lm,  2m,  and  4m 

Bottom  Volume  sound  speed  fluctuation 

3  sets 

5m/s,  15m/s,  and  45m/s 

Water  Volume  sound  speed  fluctuation 

4  sets 

0.5m/s,  lin/s,  2.5m/s,  and  5  m/s 

Table  1  MMPE  Input  Parameters 
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IV.  DATA  ANALYSIS  AND  RESULTS 


A.  MMPE  MODEL  RESULTS 

1.  Time  Domain  Reverberation  Analysis 

In  the  previous  chapter,  the  method  of  time  domain  analysis  of  the  reverberation 
loss  for  the  interfaces  and  the  bottom  volume  was  described.  To  calculate  the  structure  of 
the  reverberation  loss,  RL  from  the  interfaces  in  a  more  standard  dB  scale,  the  output 
from  the  MMPE  Model  and  post-processing  is  adjusted  by 

M't(O  =  -201og[>_,(O|],  (») 

where  p_b(t„ )  is  the  interface  (either  one)  reverberation  pressure  at  the  receiver  calculated 

in  Eq.  (94).  Similarly,  reverberation  loss  structure  from  the  bottom  volume  is  determined 
by 

RL  \,{tn)  =  -201og[|/?_v(i„)|] ,  (100) 

where  p_v(tn)  is  the  bottom  volume  reverberation  pressure  at  the  receiver  calculated  in 

Eq.  (98).  In  either  case,  the  result  is  the  predicted  reverberation  loss  structure  from  either 
of  the  interfaces  or  the  bottom  volume.  It  should  be  noted  that  this  method  predicts 
reverberation  loss  structure  and  not  actual  levels  since  we  are  neglecting  the  actual 
scattering  strengths. 

Typical  predicted  reverberation  loss  structures  for  both  interfaces  and  the  bottom 
volume  are  depicted  in  Figure  1 1 .  The  source  depth  used  for  all  of  the  numerical  analysis 
is  48m.  A  receiver  depth  of  40m  is  chosen  to  display  typical  predicted  reverberation  loss 
structures  in  this  figure  and  all  of  the  reverberation  loss  structure  plots  that  follow. 
Unless  otherwise  stated,  the  applied  source  spectrum  is  the  broadband,  SUS  spectrum 
source  with  a  center  frequency  of  350  Hz  and  a  bandwidth  of  500  Hz  calculated  over 
1024  frequency  bins. 
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Bottom  Volume.  Source  «8m.  Recerver40m 


Figure  11.  500  Hz  bandwidth,  350  Hz  center  frequency,  SUS  spectrum  reverberation  loss 

curves  for  “reference”  environment  for  source  at  48  m  and  receiver  at  40  m 

The  purpose  of  performing  this  numerical  analysis  is  to  generate  a  prediction  of 
the  expected  spatial  coherence  in  the  vertical  direction  for  an  environment  and  problem 
geometry  similar  to  that  observed  in  ASIAEX.  To  determine  the  spatial  coherence  from 
the  MMPE  Model  outputs,  a  cross  correlation  method  is  used.  In  this  method,  the 
nonnalized  cross  correlation  of  the  interface  reverberation  pressures,  p_b(t„) ,  for  the 
source  and  each  receiver  is  performed.  To  generate  a  peak  value  in  dB,  the  maximum  of 
20  log  ( | xco/r | ) ,  where  xcorr  is  the  output  of  the  normalized  cross  correlation,  is 

determined.  This  is  repeated  for  all  16  elements  and  plotted  for  both  interfaces  and  the 
volume  to  generate  an  expected  spatial  coherence. 

It  should  be  noted  that  all  correlation  curves,  unless  otherwise  stated,  are  based  on 
a  single  environmental  realization.  Ideally,  multiple  random  realizations  should  be 
computed  to  create  average  correlation  curves.  However,  due  to  the  calculation  times  of 
the  model,  such  an  approach  is  prohibitive. 

2.  Influence  of  Interface  Roughness  Variation 

Plots  of  the  predicted  reverberation  loss  structures  for  all  3  variations  in  interface 
roughness  are  plotted  in  Figure  12.  Variation  in  interface  roughness  appears  to  have  no 
strong  effect  on  the  reverberation  loss  structure  for  either  of  the  interfaces  or  the  bottom 
volume.  This  conclusion  is  consistent  with  the  results  of  previous  thesis  work^  for  the 
water/bottom  interface  and  the  bottom  volume.  A  slight  increase  in  reverberation 
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pressure  level  with  an  increase  in  the  rms  roughness  value,  which  is  particularly  apparent 
for  the  bottom/sub-bottom  interface  in  the  previous  work,  is  not  evident  here. 
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Figure  12.  Reverberation  loss  curves  for  three  variations  in  interface  roughness,  (top 
left)  water/bottom  interface;  (top  right)  bottom/sub-bottom  interface;  (bottom)  volume 


The  corresponding  peak  vertical  correlation  plots  are  depicted  in  Figure  13.  For 
both  the  water/bottom  interface  and  the  bottom  volume,  variations  in  interface  roughness 
have  very  little  effect.  The  bottom/sub-bottom  interface  displays  some  difference  in  peak 
vertical  correlation  at  the  upper  and  lower  edges  of  the  array.  A  clear  overall  trend  with 
variations  in  bottom  roughness  is  not  apparent  for  either  of  the  interfaces  or  the  bottom 
volume.  As  with  the  reverberation  loss  structure,  the  results  for  the  water/bottom 
interface  and  the  bottom  volume  are  similar  to  those  of  the  previous  thesis.  However,  a 
clear  decrease  in  vertical  correlation  with  increasing  roughness  was  previously  scen17^  for 
the  bottom/sub-bottom  interface,  and  that  trend  is  not  seen  here. 
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Figure  13.  Peak  vertical  correlation  curves  of  reverberation  loss  for  three  variations  in 
interface  roughness,  (top  left)  water/bottom  interface;  (top  right)  bottom/sub-bottom 

interface;  (bottom)  volume 


3.  Influence  of  Bottom  Sound  Speed  Perturbations 

Plots  of  the  predicted  reverberation  loss  structures  for  all  3  variations  in  bottom 
sound  speed  perturbations  are  plotted  in  Figure  14.  For  both  interfaces  and  the  bottom 
volume,  a  clear  increase  in  reverberation  loss  with  an  increase  in  bottom  sound  speed  rms 
perturbation  value  is  shown.  The  previous  thesis'-71  displayed  the  same  trend  in  all  three 
cases,  though  to  a  larger  degree.  This  trend  is  probably  caused  by  an  increase  in  the 
scatter  of  acoustic  energy  to  deeper  depths  at  shorter  ranges  caused  by  the  larger  bottom 
volume  perturbations. 
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BottonVSub- Bottom  Interface.  Source  48m  Recover  40m 
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Figure  14.  Reverberation  loss  curves  for  three  variations  in  bottom  volume  sound  speed 
perturbation,  (top  left)  water/bottom  interface;  (top  right)  bottom/sub-bottom  interface; 

(bottom)  volume 


The  corresponding  peak  vertical  correlation  plots  are  depicted  in  Figure  15. 
Except  for  the  largest  bottom  sound  speed  rms  perturbation  value,  there  is  little  change  in 
the  peak  vertical  correlation.  For  the  bottom  volume,  even  the  largest  rms  perturbation 
value  does  not  have  a  significant  effect.  Once  again,  these  results  are  similar  to  those  of 
the  previous  thesis.  The  effect  of  the  highest  perturbation  can  be  considered  an  anomaly 
since  it  produces  unrealistic  regions  where  bottom  volume  sound  speed  drops  below  that 
of  the  water  volume,  allowing  a  large  amount  of  acoustic  energy  to  transmit  directly 
across  the  interface  into  the  bottom,  where  it  is  mostly  lost.  Except  for  this  unrealistic 
case,  it  appears  that  variations  in  bottom  sound  speed  perturbations  have  little  influence 
on  peak  vertical  correlation. 
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Figure  15.  Peak  vertical  correlation  curves  of  reverberation  loss  for  three  variations  in 
bottom  volume  sound  speed  perturbation,  (top  left)  water/bottom  interface;  (top  right) 
bottom/sub-bottom  interface;  (bottom)  volume 

4.  Influence  of  Water  Volume  Sound  Speed  Turbulence 

Plots  of  the  predicted  reverberation  loss  structures  for  all  5  variations  in  water 
sound  speed  turbulence  are  plotted  in  Figure  16.  In  this  case,  only  the  reverberation  loss 
for  the  two  interfaces  was  processed.  Plots  of  the  difference  between  the  reverberation 
loss  structure  for  the  reference  environment,  which  contained  no  water  volume  sound 
speed  turbulence,  and  the  4  increased  rms  turbulence  values  are  included  in  this  figure  to 
better  represent  the  small  influence  of  increased  sound  speed  turbulence.  A  close  look  at 
these  difference  plots  shows  that  the  generally  small  difference  between  the  reference  and 
more  turbulent  environments  also  maintains  its  own  structure  over  different  turbulence 
values.  This  is  probably  an  artifact  caused  by  the  use  of  the  same  water  volume  sound 
speed  realization  modified  by  a  scaling  parameter  to  generate  the  different  turbulent 
fields. 
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Figure  16.  Reverberation  loss  curves  and  difference  plots  for  five  variations  in  water 
volume  sound  speed  turbulence,  (top)  water/bottom  interface  RL  plot  and  difference 
plot;  (bottom)  bottom/sub-bottom  interface  RL  plot  and  difference  plot 


The  corresponding  peak  vertical  correlation  plots  are  depicted  in  Figure  17.  As 
with  the  reverberation  loss  curves,  there  is  very  little  difference  in  peak  vertical 
correlation  induced  by  variations  in  water  volume  turbulence  for  either  interface.  In  fact, 
the  difference  between  variations  in  the  water/bottom  interface  is  only  a  small  fraction  of 
a  dB.  Though  the  influence  of  turbulence  was  not  accounted  for  in  the  previous  thesis, [7] 
these  results  are  very  similar  to  those  for  moderately  small  variations  in  the  SVP.  With 
this  result  in  mind,  it  is  apparent  that  changes  in  the  water  volume  sound  speed  of  this 
scale  have  little  influence  on  either  the  reverberation  loss  structure  or  the  peak  vertical 
correlation. 
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Two-Way  Rcvetberabon  Peak  Temporal  Correlation  idBi  For  Wafer  iBoCom  M  erf  ace 


Peak  Correlation  «JB) 


Peak  Correlation  (dB) 


Figure  17.  Peak  vertical  correlation  curves  and  difference  plots  for  five  variations  in 
water  volume  sound  speed  turbulence,  (top)  water/bottom  interface  RL  plot  and 
difference  plot;  (bottom)  bottom/sub-bottom  interface  RL  plot  and  difference  plot 


5.  Influence  of  Bandwidth 

Our  goal  is  to  generate  a  useful  model  for  use  as  a  basis  for  comparison  for  spatial 
coherence  analysis  of  broadband  data.  To  accomplish  this,  the  numerical  analysis  must 
include  the  peak  vertical  correlation  structure  for  both  interfaces  and  the  bottom  volume 
over  several  different  bandwidths.  In  this  case,  a  bandwidth  reduction  to  250  Hz  over 
512  frequency  bins  is  generated  directly  by  running  the  MMPE  Model  for  the  new 
spectrum.  However,  results  for  even  smaller  bandwidths  are  generated  by  removing  the 
unwanted  components  from  the  original  interface  and  bottom  volume  reverberation 
pressures, p_bv(t„),  which  generates  the  reverberation  pressures  for  these  new,  smaller 

bandwidths.  Reverberation  loss  and  peak  vertical  correlation  are  then  calculated  from 
these  reverberation  pressures  in  the  same  manner  as  before. 

Reducing  the  bandwidth  from  the  500  Hz  bandwidth,  350  Hz  center  frequency 
SUS  spectrum  in  the  reference  environment  results  in  the  peak  vertical  correlations 
depicted  in  Figure  18.  For  these  new  frequency  bands,  the  center  frequency  remained 
fixed  at  350  Hz,  but  the  bandwidth  was  reduced  to  250  Hz,  125  Hz,  and  62.5  Hz.  These 
combinations  of  bandwidth  and  center  frequency  result  in  increasing  fc/BW  ratios  of 
0.7,  1.4,  2.8,  and  5.6,  respectively.  As  opposed  to  the  example  given  in  the  correlation 
method  theory,  no  exploitable  trend  in  bottom  or  interface  vertical  peak  correlation 
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presents  itself.  In  fact,  the  only  apparent  trend  is  that  both  interfaces  and  the  volume 
display  somewhat  higher  levels  of  peak  vertical  correlation  with  decreasing  bandwidth. 


Two-W»y  Reverberation  Peak  Temporal  Correlation  (dB)  For  500  Hz  Bandwidth 


Peak  Correlation  (dB) 

Two-Way  Reverberation  Peak  Temporal  Correlation  (dB)  For  125  Mi  Bandwidth 


Peak  Correlation  (dB) 


Two-Way  Reverberation  Peak  Temporal  Correia  Son  (dB)  For  250  Hi  Bandwidth 


Peak  Correlation  (dB) 

Two-Way  Reverberation  Peak  Temporal  Correlation  (dB)  For  62  5  Hi  Bandwidth 


Peak  Correlation  (dB) 


Figure  18.  350  Hz  center  frequency,  SUS  spectrum  peak  vertical  correlation  curves  for 

four  bandwidths.  (top  left)  500  Hz  bw;  (top  right)  250  Hz  bw;  (bottom  left)  125  Hz  bw; 

(bottom  right)  62.5  Hz  bw 


6.  Influence  of  Source  Spectrum 

In  the  previous  thesis,1 ^  similar  numerical  analysis  was  conducted  using  a 
different  source  spectrum.  In  particular,  a  250  Hz  bandwidth,  250  Hz  center  frequency 
Hanning  source  spectrum  was  used.  For  this  source  spectrum,  the  peak  vertical 
correlation  decreased  much  more  rapidly  for  the  volume  than  either  interface.  This 
suggests  that  an  exploitable  difference  in  the  trend  in  bottom  or  interface  vertical  peak 
correlation  with  decreasing  bandwidth  may  be  found  in  the  numerical  analysis  if  a 
different  source  spectrum  is  used. 

The  first  component  of  the  source  spectrum  that  is  changed  to  test  this  notion  is 
the  applied  amplitude  spectrum.  A  Hanning  window  replaces  the  SUS  spectrum  and  the 
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MMPE  model  is  run  for  the  500  Hz  bandwidth,  350  Hz  center  frequency  band  in  the 
reference  environment.  All  smaller  bandwidth  calculations  are  made  in  a  similar  fashion 
as  before.  The  results  of  this  change  in  the  source  are  shown  in  Figure  19.  For  this  case, 
a  trend  with  decreasing  bandwidth  can  be  seen.  Despite  the  fact  that  peak  vertical 
correlation  doesn’t  appear  to  change  significantly  for  either  interface,  peak  vertical 
correlation  for  the  bottom  volume  is  clearly  less  than  that  of  the  interfaces  for  the  500  Hz 
bandwidth.  As  the  bandwidth  is  reduced,  peak  vertical  correlation  for  the  bottom  volume 
increases  to  the  point  where  it  appears  similar  to  the  interfaces.  This  is  an  exploitable 
trend. 


Peak  Correlation  (dB) 

Two-Way  Reverberation  Peak  Temporal  Correlation  (dB)  For  125  He  Banchrtdtti 


Peak  Correlation  (dB) 


Two-Way  Reverberation  Peak  Temporal  Coneiebon  (dB)  For  250  He  Bandwtcth 


Peak  Correlation  (dB) 

Two-Way  Reverberation  Peak  Temporal  Correlation  (dB)  For  62  5  He  Bandwidth 


Peak  Correlation  (dB) 


Figure  19.  350  Hz  center  frequency,  Hanning  spectrum  peak  vertical  correlation  curves 

for  four  bandwidths.  (top  left)  500  Hz  bw;  (top  right)  250  Hz  bw;  (bottom  left)  125  Hz 

bw;  (bottom  right)  62.5  Hz  bw 


The  most  likely  explanation  for  this  change  in  peak  vertical  correlation  with  the 
source  spectrum  lies  in  the  structure  of  the  SUS  spectrum  itself.  The  SUS  spectrum 
applied  to  the  500  Hz  bandwidth,  350  Hz  center  frequency  case  is  shown  in  Figure  20. 
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Note  that  several  narrow  peaks  rise  above  the  rest  of  the  frequency  weighting  values. 
These  peaks  are  effectively  small  bandpass  filters.  In  light  of  this,  the  mild  effect  on  peak 
correlation  observed  for  reduced  bandwidth  when  the  explosive  SUS  spectrum  is  applied 
is  understandable.  The  output  for  the  larger  bandwidths  is  essentially  passed  through 
small  bandpass  filters  already,  and  further  reduction  in  bandwidth  should  yield  little 
change. 
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Figure  20.  SUS  Spectrum  for  500  Hz  Bandwidth,  350  Hz  Center  Frequency 

7.  Influence  of  Center  Frequency 

Another  difference  between  the  source  used  in  this  thesis  and  the  previous  thesis 
is  the  center  frequency.  To  test  the  influence  of  this  difference  as  well,  the  MMPE  model 
is  run  for  a  250  Hz  bandwidth,  250  Hz  center  frequency  SUS  spectrum  band  in  the 
reference  environment.  The  125  and  62.5  Hz  bandwidth  calculations  are  then  performed 
in  the  same  manner  as  before.  A  500  Hz  bandwidth  is  not  calculated  in  this  case  since  it 
includes  the  DC  (i.e.,  0  Hz)  and  very  low  frequency  (VLF)  components  that  are 
problematic  for  the  MMPE  Model.  Note  that  this  change  in  center  frequency  results  in 
new  fc/BW  ratios  of  1,  2,  and  4,  respectively.  The  results  for  these  bands  are  shown  in 
Figure  21.  As  with  the  350  Hz  centered  SUS  spectrum,  no  discernable  trends  in  peak 
vertical  correlation  is  noted. 
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Two-Way  Reverberation  Peak  Temporal  Correlation  <d8)  For  250  Hz  Bandwidth 
Soiree  Depth  40m 


Two-Way  Reverberation  Peak  Temporal  Correlation  (d8)  For  125  Hz  Bandwidth 
Soiree  Depth  40m 


Two-Way  Reverberation  Peak  Temporal  Correlation  (dB)  Far  62  5  Hz  Bandwidth 


Figure  21.  250  Hz  center  frequency,  SUS  spectrum  peak  vertical  correlation  curves  for 

four  bandwidths.  (top  left)  250  Hz  bw;  (top  right)  125  Hz  bw;  (bottom)  62.5  Hz  bw 


The  final  source  spectrum  to  be  considered  is  the  exact  spectrum  used  in  the 
previous  thesis,  a  250  Hz  bandwidth,  250  Hz  center  frequency  Hanning  window  spectrum 
frequency  band  in  the  reference  environment.  The  results  for  this  band  are  shown  in 
Figure  22.  As  with  the  350  Hz  centered  Hanning  window  spectrum  results,  the  same 
clear  and  exploitable  trend  is  shown.  Peak  vertical  correlation  values  for  both  interfaces 
again  remain  essentially  unchanged,  while  the  peak  vertical  correlation  for  the  bottom 
volume  is  much  lower  at  larger  bandwidths  and  increases  to  levels  close  to  that  of  the 
interfaces  with  decreasing  bandwidth. 

It  is  important  to  note  that,  given  the  right  source  spectrum,  there  does  appear  to 
be  a  trend  in  the  vertical  coherence  results  which  distinguishes  interface  and  volume 
reverberation  contributions  at  large  bandwidth.  This  is  consistent  with  the  theoretical 
results  of  Ivakin.[13]  However,  in  contrast  to  the  theory,  the  numerical  analysis  suggests 
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that  a  decrease  bandwidth  tends  to  increase  the  coherence  of  the  volume  reverberation. 
The  cause  of  this  discrepancy  is  unknown,  but  possible  influences  from  the  propagation 
could  be  the  multipath  structure  and/or  the  background  sound  speed  gradient  in  the 
bottom  volume.  Analysis  of  the  impact  of  these  phenomena  remains  for  future  work. 


Two-Way  Reverberation  Peak  Temporal  Correlation  <dB)  For  250  Mr  Bandwidth  Two-Way  Reverberation  Peak  Temporal  Comeiabon  <dB)  For  125  Hz  Bandwidth 


Figure  22.  250  Hz  center  frequency,  Hanning  spectrum  peak  vertical  correlation  curves 

for  four  bandwidths.  (top  left)  250  Hz  bw;  (top  right)  125  Hz  bw;  (bottom)  62.5  Hz  bw 

B.  ASIAEX  DATA  RESULTS 
1.  General  Results 

Among  the  experiments  that  were  conducted  during  the  ASIAEX,  one  set  of 
experiments  involved  deploying  explosive  charges  from  the  research  vessel  Shiyan-3. 
The  reverberation  from  the  shallow  water  environment  was  then  recorded  by  a  vertical 
array  that  was  suspended  from  the  ship.  The  recorded  signals  form  a  data  set  from  nine 
38g  charges  and  four  1kg  charges.  Received  signals  were  digitally  recorded  at  a 
sampling  frequency  of  6000Hz. 
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To  record  the  reverberant  signals,  a  32  element  vertical  array,  designated  VLA1, 
was  employed.  This  array  spanned  the  majority  of  the  water  column,  covering  about  86m 
of  its  100m  depth.  The  spacing  of  the  elements  of  this  array  was  semi-irregular  as  shown 
in  Table  2.  Two  of  these  elements  malfunctioned  during  the  experiment,  I- 11  and  1-15, 
and  the  recordings  from  these  elements  were  not  used  in  the  data  analysis.  An  index 
number  was  assigned  to  the  remaining  30  working  elements,  which  is  used  in  some  of  the 
figures  that  follow. 


Hydrophone 

Depth 

Spacing 

Remarks 

Hydrophone 

Depth 

Spacing 

Remarks 

1-32 

4.60m 

- 

Index  30 

1-16 

32.50m 

2.00m 

Index  14 

1-31 

6.33m 

1.73m 

Index  29 

1-15 

36.50m 

4.00m 

Bad  element 

1-30 

8.13m 

1.80m 

Index  28 

1-14 

40.50m 

4.00m 

Index  13 

1-29 

10.21m 

2.08m 

Index  27 

1-13 

44.50m 

4.00m 

Index  12 

1-28 

11.69m 

1.46m 

Index  26 

1-12 

48.50m 

4.00m 

Index  11 

1-27 

13.35m 

1.66m 

Index  25 

1-11 

52.50m 

4.00m 

Bad  element 

1-26 

15.02m 

1.67m 

Index  24 

1-10 

56.50m 

4.00m 

Index  10 

1-25 

16.81m 

1.79m 

Index  23 

1-9 

60.50m 

4.00m 

Index  9 

1-24 

18.61m 

1.80m 

Index  22 

1-8 

64.50m 

4.00m 

Index  8 

1-23 

20.38m 

1.77m 

Index  2 1 

1-7 

68.50m 

4.00m 

Index  7 

1-22 

21.96m 

1.58m 

Index  20 

1-6 

72.50m 

4.00m 

Index  6 

1-21 

23.68m 

1.72m 

Index  19 

1-5 

76.50m 

4.00m 

Index  5 

1-20 

25.36m 

1.68m 

Index  18 

1-4 

80.50m 

4.00m 

Index  4 

1-19 

27.06m 

1.70m 

Index  17 

1-3 

82.50m 

2.00m 

Index  3 

1-18 

28.75m 

1.69m 

Index  16 

1-2 

86.50m 

4.00m 

Index  2 

1-17 

30.50m 

1.75m 

Index  15 

1-1 

90.50m 

4.00m 

Index  1 

Table  2  Element  Spacing  and  Depth  for  Vertical  Array  VLA1  Used  in  ASIAEX 

An  unexpected  difference  in  the  recorded  signals  between  different  elements  was 
discovered  early  in  the  analysis.  The  recorded  signals  from  the  lower  16  elements  of 
VLA1,  1-1  through  1-16,  displayed  a  different  received  signal  shape  than  signals  from  the 
upper  16  elements,  1-17  through  1-32.  An  example  of  this  difference  from  the  first  38g 
charge  for  the  two  adjacent  elements  of  these  two  16  element  sets  is  shown  in  Figure  23. 
Further  analysis  shows  that  the  signals  from  the  upper  elements  contain  stronger  low 
frequency  components  than  the  signals  recorded  by  the  lower  elements.  Element 
placement  within  the  water  column  should  not  account  for  this  difference,  as  1-16  and  I- 
17  are  only  2m  apart.  Instead,  it  has  been  confirmed  that  the  settings  for  these  two  16 
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element  sets  were  indeed  different,  though  the  exact  nature  of  this  difference  remains 
unknown. 


Time  <*ec|  Time  (seel 

Figure  23.  Two  Received  Signals  From  the  First  38g  Charge. 


Plots  for  the  received  signals  recorded  by  the  same  two  elements,  1-16  and  1-17, 
for  the  first  1kg  charge  are  shown  in  Figure  24.  The  difference  in  the  shape  of  the 
received  signals  from  the  two  halves  of  the  array  is  still  quite  evident.  These  plots  also 
reveal  that  the  signals  recorded  from  this  charge  display  a  large  amount  of  clipping.  This 
was  present  in  all  four  of  the  1kg  charge  recordings.  Some  signal  clipping  was  also 
evident  in  the  recorded  signals  for  the  final  three  38g  charges. 


Received  Signal  at  Element  M6 


Received  Signal  at  Element  1-17 


Figure  24.  Two  Received  Signals  From  the  First  1kg  Charge. 


Vertical  peak  correlation  for  the  recorded  data  is  determined  in  a  manner  similar 
to  that  used  for  the  numerical  analysis  results.  In  this  case  there  is  no  clear  reference 
element  since  none  of  the  elements  corresponds  directly  to  the  source.  Instead,  the 
normalized  cross  correlation  of  all  30  working  elements  is  performed  with  respect  to  each 
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of  the  others.  The  maximum  value  of  each  of  these  cross  correlations  is  determined  to 
form  a  30x30  matrix  of  peak  correlation  values.  These  peak  normalized  cross  correlation 
values  are  then  converted  to  a  dB  scale  in  the  same  manner  as  before,  201og(|xco/r|) . 

Two  methods  are  used  to  display  the  peak  vertical  correlation  results  for  the 
ASIAEX  data.  Due  to  the  semi-irregular  spacing  of  VLA1,  the  relative  spacing  between 
elements  rarely  aligned,  making  a  spatial  average  impossible  to  determine.  Instead,  the 
plots  of  peak  vertical  correlation  with  respect  to  all  30  working  elements  are 
superimposed  on  the  same  plot  against  an  axis  of  the  relative  depth  difference  between 
elements  for  each  correlation  calculation.  To  visually  sort  this  wealth  of  information,  an 
image  plot  is  also  generated  where  the  peak  correlation  matrix  is  displayed  as  an  image 
with  color-coding  representing  the  numerical  correlation  value  in  dB.  Note  that  the 
autocorrelation  values  run  diagonally  through  the  image.  The  axes  of  this  plot  use  the 
previously  defined  element  index  number  from  Table  2.  An  example  of  both  plot  types 
for  the  first  38g  charge  is  shown  in  Figure  25. 

Pe*fc  Vert*c»l  Cross  Correlation  (dB  i  Ful  3000  Hz  BanrtiwdP)  Pss*  Vertical  Cross  CorreMwv  Fui  3000  Hz 


Figure  25.  Peak  Vertical  Correlation  Between  Elements  For  38g  Charge 

The  effect  of  the  difference  in  the  recorded  signals  between  the  upper  and  lower 

portions  of  the  array  can  be  clearly  seen  in  the  image  plot  in  Figure  25.  Peak  vertical 

correlation  values  are  much  higher  between  elements  of  the  same  half  of  the  array  than 

they  are  between  elements  in  different  halves  of  the  array.  In  further  analysis,  the 

received  signals  are  filtered  using  a  bandpass  filter  prior  to  conducting  the  cross 

correlations.  For  bands  that  contain  DC  and  the  low  frequency  components,  the  bi-polar 

structure  seen  in  the  full  spectrum  analysis  is  still  present.  For  bands  that  remove  these 
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low  frequency  components,  the  bi-polar  structure  is  absent  and  peak  vertical  correlation 
values  are  generally  lower  overall  (though  significantly  higher  between  elements  located 
in  different  halves  of  the  array).  Examples  of  this  trend  are  evident  when  Figure  25  and 
26  are  compared. 


Peak  Vertical  CorreJatton  0-1 000 Hi 


Figure  26.  Peak  Vertical  Correlation  Between  Elements  For  38g  Charge,  0-1000  Hz 

Band 


Pert  Vertical  Correlation  ufi)  500-1500HT 
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Peak  Vertical  Correlation  500- 1500 Hi 


Figure  27.  Peak  Vertical  Correlation  Between  Elements  For  38g  Charge,  500-1500  Hz 

Band 


To  this  point,  only  examples  for  the  38g  charges  have  been  displayed.  The  same 

plot  types  as  shown  above  are  displayed  for  the  first  1kg  charge  in  Figures  28,  29,  and  30. 

In  general,  the  same  trends  are  apparent  for  the  larger  explosive  charges  as  were  seen  for 

the  38g  charges.  Overall  peak  vertical  correlation  values  are  much  lower  for  these  1kg 

charges  than  for  their  38g  counterparts.  This  is  probably  a  result,  at  least  in  part,  of  the 

large  amount  of  signal  clipping  observed.  The  high  correlation  between  elements  in  the 

upper  portion  of  the  array  relative  to  the  elements  of  the  lower  portion  of  the  array  for 

bands  that  contain  DC  and  low  frequencies  is  generated  by  the  differences  in  the  recorded 
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signals  in  these  two  halves  of  VLA1.  Due  to  this  low  correlation,  the  remaining  analysis 
will  focus  on  the  38g  charges. 


Peak  Vertical  Croat  Correction  (dBl  Pul  3000  HU  Bandcdm 


Peak  Vertical  Cross  Correction  FVI 3000  HU  BandMdth 


Figure  28.  Peak  Vertical  Correlation  Between  Elements  For  1kg  Charge 


Peak  Vertical  Correction  (dBl.  CM 000 HU 


Peak  Vertical  Correction  0-1000H* 


Figure  29.  Peak  Vertical  Correlation  Between  Elements  For  1kg  Charge,  0-1000  Hz 

Band 


Peek  Vertical  Correlation  500- 1500 HU 


Figure  30.  Peak  Vertical  Correlation  Between  Elements  For  1kg  Charge,  500-1500  Hz 

Band 
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A  closer  analysis  of  the  38g  charges  reveals  that  each  charge  displays  particularly 
high  peak  vertical  correlation  values  in  the  vicinity  of  100  Hz.  For  each  charge,  the 
specific  center  frequency  for  which  a  band  displays  the  highest  peak  vertical  correlation 
is  given  in  Table  3.  No  strong  environmental  source  is  observed  at  or  near  these 
frequencies  prior  to  the  recorded  explosion.  The  cause  of  this  region  of  high  vertical 
correlation  could  be  an  artifact  of  the  array,  the  recording  equipment,  or  an  oddity  in  the 
environment.  An  attempt  to  remove  the  explosive  source  spectrum  from  the  recorded 
signal  using  a  similar  SUS  spectrum  generated  in  the  same  manner  as  in  the  numerical 
analysis  failed  to  remove  this  region  of  high  vertical  correlation,  but  this  result  is  far  from 
conclusive.  It  should  be  noted  that  the  numerical  analysis  does  not  include  frequencies  of 
100  Hz  or  below,  so  the  effect  of  this  high  correlation  region  on  comparison  between  data 
and  numerical  analyses  is  minimal. 


Charge  Number 

High  Correlation 
Center  Frequency 

Charge  Number 

High  Correlation 
Center  Frequency 

1 

80  Hz 

6 

90  Hz 

2 

105  Hz 

7 

85  Hz 

3 

105  Hz 

8 

90  Hz 

4 

80  Hz 

9 

95  Hz 

5 

80  Hz 

Table  3  High  Peak  Vertical  Correlation  Center  Frequencies  for  38g  Charges 

2.  Influence  of  Bandwidth  on  Peak  Vertical  Correlation 

The  results  of  the  numerical  analysis  showed  that,  for  a  broadband  explosive 
(SUS)  source,  there  was  no  exploitable  trend  with  decreasing  bandwidth  observed 
between  the  peak  vertical  correlation  values  from  reverberation  caused  by  the  interfaces 
and  that  caused  by  the  bottom  volume.  In  fact,  the  only  trend  in  peak  vertical  correlation 
with  decreasing  bandwidth  is  that  overall  correlation  values  become  somewhat  higher  as 
bandwidth  is  decreased.  This  trend  is  clear  in  the  350  Hz  centered  band  results. 

To  compare  the  results  of  the  numerical  analysis  to  the  recorded  data,  peak 
vertical  correlation  calculations  are  detennined  for  the  precise  frequency  bands  used  in 
the  numerical  analysis.  In  each  case,  bandwidth  is  varied  in  an  attempt  to  verify  the  trend 
of  increasing  correlation  with  decreasing  bandwidth  indicated  by  the  numerical  analysis. 
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The  first  of  these  bands  to  be  reproduced  is  the  set  centered  at  350  Hz.  As  in  the 
numerical  analysis,  bandwidths  of  500  Hz,  250  Hz,  125  Hz,  and  62.5  Hz  are  used.  The 
results  of  this  analysis  for  the  first  38g  charge  are  displayed  in  Figures  31  and  32.  As 
with  the  two  previous  sets  of  smaller  frequency  bands,  this  analysis  supports  the 
predicted  trend  of  increasing  correlation  with  decreasing  bandwidth. 


Peak  Vertical  Conett&on  <d8).  500  Hi  BW  centered  on  350  H* 


Peak  Vertical  ConettWn  idBi.  250  Hj  BW  centered  on  350  Hj 


Element  Index 


Figure  3 1 .  Peak  Vertical  Correlation  Image  Plots  For  38g  Charge,  350  Hz  Center 

Frequency  Bands 
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P*«k  Vertical  Correlation  500  Hi  BW  centered  on  350  Hi 


Peak  Vertical  Correlation  250  Hi  BW  centered  on  350  Hi 


Relative  Depth  (ml 


Peak  Vertical  Correlation.  62  5  Hi  BW  centered  on  350  Hi 


Relative  Depth  (m| 


Figure  32.  Peak  Vertical  Correlation  Plots  For  38g  Charge,  350  Hz  Center  Frequency 

Bands 


The  final  set  of  frequency  bands  to  be  considered  are  those  centered  at  250  Hz. 
As  in  the  numerical  analysis,  bandwidths  of  250  Hz,  125  Hz,  and  62.5  Hz  are  used.  The 
results  of  this  analysis  for  the  first  38g  charge  are  displayed  in  Figures  33  and  34.  Once 
again,  this  analysis  supports  the  predicted  trend  of  increasing  correlation  with  decreasing 
bandwidth,  as  does  all  of  the  data  analysis  conducted.  Because  this  trend  was 
consistently  more  noticeable  in  the  numerical  analysis  for  the  volume  reverberation,  the 
result  of  this  data  analysis  may  indicate  a  dominant  contribution  from  the  bottom  volume 
fluctuations. 
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Peak  Verfccel  Correlation  (dBi.  250  Hi  BW  certered  on  250  Hi 


Peak  Verted Correlation  <dBi  125  Hi  BW  certered  on  250  Hz 


Berner*  Index  Berner*  Index 


•1-10 


H-12 


Figure  33.  Peak  Vertical  Correlation  Image  Plots  For  38g  Charge,  250  Hz  Center 

Frequency  Bands 


62 


8 

|  -8 


-16 

•16 


•^oo 


Figure  34. 


Peak  Vertical  Correlation.  250  Hi  BW  centered  on  250  m 
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Peak  Vertical  Correlation.  62  5  Hi  BW  centered  on  250  Hi 


Peak  Vertical  Correlation  Plots  For  38g  Charge,  250  Hz  Center  Frequency 

Bands 
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V.  SUMMARY 


The  focus  of  this  thesis  was  to  examine  the  spatial  coherence  of  broadband 
reverberation  data  in  a  shallow  water  environment  using  explosive  charges,  and 
determine  if  unique  features  of  the  interface  or  volume  reverberation  structure  could  be 
used  to  distinguish  the  dominant  scattering  mechanism.  To  accomplish  this,  a  numerical 
model  was  employed  to  predict  the  peak  vertical  correlation  structure  for  the 
water/bottom  interface,  the  bottom/sub-bottom  interface,  and  the  bottom  volume  over 
several  different  bandwidths  in  this  environment.  The  environment  and  geometry  were 
chosen  to  approximate  the  conditions  present  during  the  East  China  Sea  component  of 
ASIAEX  and  to  match  the  conditions  used  in  previous  work  as  much  as  practical. 
Included  in  this  numerical  analysis  was  the  incorporation  of  an  explosive  SUS  charge 
frequency  spectrum,  a  change  from  the  Hanning  window  spectrum  used  in  previous 
work,  to  better  model  the  explosive  sources  used  in  ASIAEX.  A  larger  500  Hz 
bandwidth  calculated  over  1024  frequency  bins  was  also  chosen  to  more  accurately 
simulate  a  wideband  source  than  the  previous  250  Hz  bandwidth  calculated  over  512 
frequency  bins. 

Several  environmental  parameters  were  also  varied  with  respect  to  a  set  of 
reference  conditions  to  determine  their  influence  on  the  peak  vertical  correlation 
structure.  The  environmental  variations  considered  were  variations  in  the  degree  of 
roughness  for  the  water/bottom  and  the  bottom/sub-bottom  interfaces,  variations  in  the 
magnitude  of  bottom  volume  sound  speed  perturbations,  and  variations  in  the  magnitude 
of  water  volume  sound  speed  turbulence.  These  variations  in  the  environment  displayed 
little  overall  influence  on  the  peak  vertical  correlation  structure.  In  general,  the  few 
variations  that  were  observed  were  not  strong  and  lacked  a  clear  trend.  The  one  instance 
where  a  clear  trend  was  observed  was  an  increase  in  peak  vertical  correlation  for  the 
largest  bottom  volume  sound  speed  perturbation.  However,  this  only  occurred  when 
given  an  unrealistic  magnitude  of  the  volume  perturbation,  which  resulted  in  a  significant 
drop  in  propagation  energy,  and  should  not  be  considered  a  useful  result. 
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The  numerical  analysis  of  the  modeled  explosive  source  data  determined  that  no 
exploitable  trend  in  peak  vertical  correlation  existed  with  decreasing  bandwidth.  Peak 
correlation  values  from  both  interfaces  appeared  roughly  equal  to  those  from  the  bottom 
volume  in  all  cases.  Only  an  increasing  trend  in  the  peak  vertical  correlation  for  all  three 
of  the  contributors  was  observed.  Analysis  of  the  recorded  data  collected  during 
ASIAEX  was  consistent  with  the  results  of  this  numerical  analysis.  In  four  separate 
cases,  including  the  2  cases  specifically  modeled  in  the  numerical  analysis,  a  decrease  in 
bandwidth  resulted  in  an  increase  in  vertical  correlation.  However,  since  no  exploitable 
trend  was  observed  in  the  numerical  analysis  to  distinguish  the  reverberation  mechanism, 
no  further  information  regarding  interface  or  volume  scattering  could  be  determined  from 
the  recorded  data. 

The  previous  work  suggested  that  an  exploitable  trend  in  peak  vertical  correlation 
with  decreasing  bandwidth  does  exist  for  this  environment  and  geometry,  at  least  when  a 
different  source  spectrum  is  used.  This  was  confirmed  by  performing  the  numerical 
analysis  using  the  source  from  the  previous  work,  a  250  Hz  bandwidth,  250  Hz  center 
frequency,  Hanning  windowed  source.  The  peak  vertical  correlation  from  both  interfaces 
remained  constant  with  decreasing  bandwidth,  but  peak  vertical  correlation  from  the 
bottom  volume  increased  from  levels  clearly  below  that  of  the  interfaces  to  levels  similar 
to  those  of  the  interfaces.  This  qualifies  as  an  exploitable  trend  which  spatial  analysis  of 
recorded  signals  could  leverage  to  discriminate  interface  scatter  from  volume  scatter. 

Unfortunately,  the  cause  of  this  difference  between  a  lack  of  exploitable  trend  in 
the  current  work  and  its  presence  in  the  previous  work  was  determined  to  be  the  structure 
of  the  source  spectrum  applied.  The  various  narrow  peaks  present  in  the  frequency 
spectrum  of  an  explosive  source,  such  as  a  SUS,  act  like  several  small  bandpass  filters  in 
the  frequency  domain.  Since  the  source  spectrum  already  has  the  same  effect  as  running 
the  signal  through  several  narrow  bandpass  filters,  further  bandpass  filtration  of  the  signal 
has  little  effect  and  should  not  be  expected  to  display  a  trend.  As  a  result,  the  use  of  a 
broadband  explosive  spectrum  was  found  to  be  a  poor  choice  for  the  application  of  a 
spatial  correlation  method  to  discriminate  between  interface  and  bottom  volume 
scattering. 
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Despite  the  fact  that  it  did  not  yield  an  exploitable  trend,  the  numerical  analysis 
did  show  that  peak  vertical  correlation  structure  was  altered  little  by  both  direct  and 
indirect  environmental  influences.  This  suggests  that  exploitable  trends  in  spatial 
coherence  are  robust  in  the  face  of  environmental  variability  and  that  it  is  the  nature  of 
the  source  spectrum  that  is  the  primary  discriminating  factor  in  the  ability  to  yield  an 
exploitable  trend.  It  should  be  noted  that  in  the  example  provided  in  Ivakin’s  work,  the 
observed  trend  was  that  interface  coherence  was  lower  for  the  larger  pulse  length  while 
the  volume  coherence  remained  unchanged.  The  analogous  trend  with  decreasing 
bandwidth  seen  in  the  numerical  analysis  for  a  Hanning  windowed  source  showed  an 
opposite  trend.  The  volume  coherence  was  greater  for  smaller  bandwidths  (larger  pulse 
lengths)  while  the  interface  coherence  remained  unchanged.  However,  both  still  show 
that  the  difference  between  interface  and  volume  coherence  was  greatest  for  larger 
bandwidths  (shorter  pulse  length).  The  reason  for  the  difference  in  trend  between  theory 
and  numerical  analysis  is  unclear,  but  it  should  be  noted  that  the  results  from  the  recorded 
data  displayed  a  similar  overall  trend  to  the  combined  trends  from  the  Hanning  windowed 
source  analysis,  such  that  the  overall  vertical  coherence  increased  for  smaller  bandwidths. 
Therefore,  there  is  some  indication  that  the  volume  rather  than  the  interface  may 
dominate  reverberation  in  the  measured  data,  although  the  use  of  an  explosive  spectrum 
may  preclude  such  a  conclusion. 

With  the  conclusion  of  this  thesis,  recommendations  for  future  work  include  an 
investigation  of  the  effects  of  both  rough  sea  surface  scatter  and  multiple  radial  returns. 
The  incorporation  of  the  effects  of  rough  sea  surface  scatter  into  the  MMPE  Model  and 
the  investigation  of  its  predicted  influence  on  peak  vertical  correlation  would  add  the 
effects  of  variation  in  another  major  environmental  component  on  peak  vertical 
correlation  structure  to  our  knowledge  base.  An  investigation  of  the  influence  of  multiple 
radial  returns  on  the  signal  coherence  would  make  this  type  of  analysis  even  more 
realistic.  Further  work  could  also  be  perfonned  to  detennine  the  discrepancy  between  the 
theoretically  predicted  influence  of  bandwidth  and  the  observed  trends  in  the  numerical 
analysis. 

If  further  analysis  supports  the  conclusion  that  volume  scatter  dominates  the 

ASIAEX  reverberation  data,  then  it  may  be  possible  to  combine  measured  and  modeled 
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results  to  invert  for  the  volume  scattering  strength,  Sv .  Evaluation  of  this  could  then  lead 

to  an  estimate  of  the  volume  perturbation  spectral  strength,  which  could  be  compared  to 
results  obtained  from  bottom  core  samples  taken  from  the  region.  Furthennore,  such  a 
successful  discrimination  of  the  volume  reverberation  from  active,  broadband  returns 
could  be  exploited  in  future  processing  algorithms  to  aid  in  the  detection  of  other 
scattering  targets  of  interest. 
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